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ABSTRACT 


This work reviews the existing literature on slam force 
prediction methods and discusses some of the theory of ship 
bottom slamming. One of these prediction methods is chosen 
as the best and further improvements are made on it by uSing 
a new model of slamming. 


The work is designed to determine a slam forcing function 
Suitable for studying ship hull girder dynamic response. 
Slam forces are determined by integrating the pressure time 
history over the portion of the hull subject to slamming. 
A computer program iS included which determines the slam 
forces and impulses from the ship's lines, principal dimen- 
Sions, and seakeeping data. The method is applied to two 
ships, the MARINER and the FOTINI-L, and the results 
presented. 
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aes INTRODUCTION 

The primary objective of this thesis is to review the 
existing literature on slamming and to cevelop a method of 
predicting slam forces to determine hull girder response. 

This project iS a portion of an overall study presently under- 
way at M.I.T. sponsored by the American Bureau of Shipping 
(A.B.S.) on hull stiffness criteria. 

Presently, explicit criteria for ship stiffness do not 
exist. The ship classification societies have controlled stiff- 
ness or deflections by limiting allowable Length to Depth (L/D) 
ratios [1]*. With the introduction of new materials such as 
high strength steels and aluminum in ship construction, the 
question of stiffness and what should be allowed becomes more 
important. 

Considering the present L/D ratios as hull girder stiff- 
ness criteria, what other factors that affect stiffness could 
be used in design? Bending stiffness of a ship's hull can be 
varied within reasonable limits by the following means: 

1. By changing the basic hull depth, 

2. By deSigning to higher stress levels with higher 

Strength steels, and 

3. By uSing materials with different elastic moduli. 

Since other factors can more properly be used to control stiff- 


ness other than L/D, it is not satisfactory to use this as a 


*Numbers in square brackets deSignate References at the end 
of the paper. 
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criterion for design. Probably a new criterion based 
specifically on a measure of stiffness should be used. 

Before any quantitative information could be developed 
for limiting hull stiffness, factors which would be changed 
by stiffness variations were assessed [2]. Assuming strength 
criteria are maintained, stiffness variations will result in 
changing factors which determine hull girder dynamic response. 
In the Reference [2] report it was concluded that the factors 
of ship stiffness which should be of concern were: 

1. Slamming response for its contribution of stress 

components, especially amidships, 

2. Springing for its cumulative stress effects, and 

3. steady-state propeller-excited vibratory motions 

and their deleterious effects upon personnel and 

main machinery components. 
Although propeller-excited vibratory motions will change with 
hull stiffness variations, it is felt that these problems can 
be dealt with by more efficient and less costly methods than 
changing the hull girder stiffness. Therefore the primary 
factors of concern are slamming and springing response with 
variations in hull stiffness. 

In the overall study the Kline-Clough lumped-mass vibra- 
tion program, S.H.V.R.S., has been used to analyze the ship 
response to slamming and springing [3]. Several ships have 
been looked at in great detail. Some of the results have been 


reported in references [3] and [4]. To date, in this work on 


RZ 





slamming, a constant slam impulse of 100 Ton-seconds has been 
used in the response analysis modeled as a half sine pulse 
applied at a single location. Results, such as stress levels, 
other than trends in response due to changes in stiffness, 
could not be determined until a better model of the slam force 
and a satisfactory means of evaluating force magnitudes from 
ship to ship could be found. 

The purpose of this work is to develop a model of the 
bottom impact slam force compatible with the S.H.V.R.S. program 
which will require a minimum of manhours and computer time yet 
Gueve adequate results. The bottom impact slam force developed 
here 1s a function of hull form, ship speed, and sea stete. 

A computer program is developed which predicts the most 
probable slam force as well as an extreme force to be used for 
design. Other slam information is also available such as 
frequency of slamming and slam pressures useful in designing 


mul plating. 
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Ti. BOTTOM IMPACT SLAMMING 

For this study the bottom impact slam is defined as the 
“heavy blow" resulting from the forebody hitting the water 
surface at a high relative velocity. This is in contrast to 
the phenomenon of bow flare slamming which is a consequence 
of emersion of the bow flare and the resultant hydrodynamic 
and hydrostatic forces. Much of the work to date has been in 
the form of predicting the local pressures developed to study 
plating damage. That can be considered the micro-problem. In 
this work we are concerned with the macro-problem which in- 
volves hull vibrations and possibly large midship stresses. 

Oakley [5] has analyzed the micro-problem in great detail. 
There are many factors that have been observed to measurably 
affect the pressure history. These factors include the 
following: shape, velocity, and phasing of the surface waves; 
geometry of the hull and hull reaction characteristics; air 
and water compressibility; two and three-dimensional air and 
water flow; bubbles, spray, and even such things as marine 
growth. Due to the complexity of the problem and the random- 
ness of the sea, model testing and statistics play an important 
role in the solution to the slamming problem. 

Tick [6] developed the equation for the frequency of 
slamming which is dependent on the following: 

1. forefoot emergence statistics, 

2. relative velocity statistics, and 

3. the angle between the keel line and the wave 


surface at the instant of impact. 
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In most work to date the equation is simplified to remove the 
angle dependence and only the regquirements of (1) forefoot 
emergence and (2) a relative velocity greater than a threshold 
velocity are considered necessary for slamming. Neglecting 
this angular dependence results in overpredicting the frequency 
of slam impact as well as the slam force. A method of dealing 
with the angular dependence is presented in Section IV of this 
work. The two requirements for slamming are met by the ship 
in Figure l. Using Tick's equation without angular dependence 
allows prediction of the frequency of slamming and is the basis 
of slam calculations in seakeeping programs [7]. The one 
necessary input to seakeeping programs is the threshold velo- 
city. How this is determined and the value to be used is 
discussed in Section IV. 

From Tick's work it is possible to conclude whether or 
not a ship iS subject to slamming under a particular set of 
operating conditions, in particular the ship speed, sea state 
and draft at the station under consideration. Given that 
slamming does occur, the next step is to determine the slam 
force. Two basic methods have been used in the past: (ir) 
integration of the slam pressure time histories over the slam 
area; and (2) determine an analytical solution based on the 
rate of change of momentum. In Section III both methods are 


discussed and one selected for this study. 


Ls 





FIGURE 1 - SHIP SLAM CRITERIA 





PHASE (1) 2HASE (II) PHASE (III) 





FIGURE 2 - PHASES OF SLAMMING 
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Before comparing the two methods, some discussion is in 
order with respect to the mechanism of slamming. Oakley [5] 
considered the time history of slamming to be made up of 
three phases: 

1. Phase (I) is the time in which the body is 
approaching the free surface up until it makes 
contact. In this phase the air flow and the 
surface wave deflection are of predominant 
importance. 

2. Phase (II) is the period when the body goes 
through the free surface and involves the problem 
of "contact" and more or less complete wetting of 
the hull. 

3. Phase (III) is the fully wetted problem in which 
the rate of change of added mass and the local 
surface deflection are dominant. 

A very flat body or "U" form hull may accelerate the free 
surface significantly during phase (I) while a more "Vv" 
Shaped section may develop greatest forces only during phase 
(III). Figure 2 shows the relationship between the three 
phases. Much of the work to date has dealt with only one or 
two of these phases at a time. This is especially true with 
the many two and three-dimensional drop tests performed. 
Chuang [8] and others have demonstrated the important role 
of the air in accelerating the free surface in phase (I) and 


as a cushioning mechanism in phase (II). The wedge-entry 


ny, 





problem has been studied by many authors, most recently by 

Chuang [9], and agreement for large deadrise angles (phase 

(III) only) is good. The theory breaks down for low deadrise 

angles where phase (I) and phase (II) become important. 
Therefore any solution to the slam forcing function 

should take into account the forces developed during all 

three of the phases unless one phase can be shown to be much 

Teecieemeenan the Others. At this time it appears that neglect- 


ing phase (I) and phase :II) is only valid for "V" form sections. 
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mits CLAM FORCE PREDICTION METHODS 

During the last twenty years a great deal has been written 
on slamming, attesting t» the difficulty of a satisfactory 
solution to the problem. Much of the work has been done along 
with ship vibration anal’sis techniques. This study does not 
go into the question of structural response. It is felt that 
adequate methods exist for analyzing the structural response 
of a ship if the various forcing functions can be adequately 
represented. The S.H.V.R.S. program described in Reference [3] 
which has been used in this overall study has been shown to 
give satisfactory results. Therefore the model of slam force 
was chosen with the S.H.V.R.S. program in mind. 

The slam force prediction methods are of two types: 

1. The slam force is considered as resulting from 
the rate of change of momentum of the added 
mass of the ship with a buoyancy correction 
term. This method will be called the momentum 
method. 

2. The second method consists of determining the slam 
pressure and integrating spatially as well as 
temporally over the slam area. This method will 
be designated the slam pressure method. 

Both of these methods were considered and the second was 
considered more valid ani more easily applicable to the 


overall stiffness study. 
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RATE OF CHANGE OF MOMENTU4 METHOD 

As early as 1956 Szebehely [10] proposed that the slam 
force waS a result of the rate of change of momentum of the 
added mass in the region of the slam. Further studies by 
Leibowitz [ll] used this model of slamming to analyze the 
slam forces on a Dutch destroyer. Good results were obtained 
between the theory and full scale results for bow flare slams. 
More recently Mansour ano d'Oliviera [12], [13] have utilized 
the same slam force function to predict slam bending moments 
in regular waves. 

The slam force (P) is considered as the sum of an inertial 


term (P,) and a buoyancy term (P.) and is given by 


E(x; t } Py (x,t) + FP. (x,t) 


D 
= (mw) + cgA (1) 


3 
= (mw) + 


ax 
ee — + 
xe (mw) ogA 


alt 


Aha? 


x 


where m = m(x,t) is the edded mass per unit length, Ww. is the 
relative vertical velocity between the wave and the ship, and 
A is the sectional area. The total derivative is used (=) 
Since the rate of change of momentum occurs not only temporally 
DiwemalsOouspatially. 

It should be noted that this is the phase (III) problem 
described by Oakley. No account is made for the force which 


comes during phase (I) and phase (II) of the slam process. It 
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is expected that this method should show good results for "Vv" 
form hulls and bow flare slamming only. The major drawback 
to this approach, even if it were valid for the bottom impact 
slam case, is it requires a large amount of seakeeping infor- 
mation to determine m(x,t). In reality the added mass must 
be determined as a function of draft at each station and 
frequency of encounter fcr each wave considered. To obtain 
this information requires a large expenditure in manhours and 
computer time. The end result is also valid only for regular 


waves. 


Sim PRESSURE PREDICTION METHOD 

A great deal of the initial work in this area has been 
done to solve the micro-}; roblem--that is, to determine the 
slam pressure locally fo1 plate design information. Much of 
the work in this area ha: been done by M.D. Ochi, K.M. Ochi, 
and S.L. Chuang at the Nieval Ship Research and Development 
Center. A summary of the slamming studies to date appears in 
Reference [14]. In this work by M.K. Ochi and Motter they 
develop the basic theory behind the approach to slamming pre- 
diction in this method. 

In this approach the pressure is considered as a function 
of velocity of impact and section shape only. The slam 


pressure 1S given by 


p = kt- (2) 


Za 








where k 1S a parameter to be determined from section form 
and r is the velocity of impact. 

From the many two and three-dimenSional drop tests as 
well as full scale measur2ments and model tests it has been 
found that equation (2) holds quite well for slam impact pro- 
blems. Therefore once th: form coefficient (k) can be deter- 
mined for each section, tie slam impact pressure can be 
determined if the velocit’ of impact is known. Since relative 
velocities are easily det:2rmined from seakeeping programs, 
knowing the value of k for each section allows for the deter- 
mination of the slam pressure which can be integrated over the 
hull form to determine the slam force. 

This second approach has the advantage that the form 
coefficient k is obtained from model data which includes the 
contribution from phase (I), (II), and (III) forces. In 
addition, recent work by Loukakis and Chryssostomidis [15] 
make available the seakeeping information necessary for slam 
studies for ships that can be approximated by Series 60 hull 
forms. 

The slam pressure approach was chosen since it requires 
a minimum of seakeeping information, is applicable to all 
hull forms, and the results are applicable in irregular seas. 
Unlike most engineering problems it also turns out that this 


approach is easiest to apply. 
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iv PRINCIPLES AND PROCEDURE OF THE SLAM PRESSURE METHOD 

Ochi and Motter's work described in Reference [14] pre- 
sents the theory and assumptions of the slam pressure method. 
In trying to apply the method it appeared that some of their 
assumptions resulted in overpredicting the slam force and im- 
pulse. In this section the theory and assumptions made by 
Ochi and Motter are discussed and modified where the assump- 
tions were not felt fully justified. The following topics 
will be discussed in developing the method: 

1. Frequency of slan impact, 

2. Slam pressure ard hull form coefficient, 

3. Most probable extreme pressure, 

4. Extreme pressure for design, 

5. Sea state data end approximation for storm duration, 

6. Slam impact aree of the hull, and 


7. Slam impact force and impulse. 


FREQUENCY OF SLAM IMPACT 

Tick [6] derived the equation for the frequency of slam 
impact at the bow which is dependent on the relative velocity, 
relative motion, and angle between the keel and the wave sur- 


face. By integrating over all angles to neglect the angular 


dependence, 
Rt 
il f ) 
No = 55 ar Pr {slam impact} 
r 
R! 2 ee 
me r oe Lx 
Me on vy RT Oke (ay + RD 2) 





where N = Number of slams per second 
Tee— Shipedrart dt station under consideration 
r, = Threshold velocity 
R' = Twice variance of relative motion 
Re = Twice variance of relative velocity 
The number of impacts in a storm of duration t hours, N(t), 


is then given by 


3600 t Be 7 2 
N(t) = 7 rr exp ek aia A! (4) 
ig a i 


Different values of the threshold velocity below which 
Slamming does not occur have been used, but following Ochi's 


work [14], 


2 ee Ye Le 
my = 12/7 =55 O25 3 ey o c 


will be used where f, = 12 feet/sec for a 520 FT ship and the 
Froude scaling law is applied. Because some of the slams 
predicted by equation (4) will occur above the angle which will 
cause slamming, equation (4) will overpredict the number of 
slams. This number is used to determine the extent of slam 
area and slam pressure. Overestimating N(t) will yield a 
larger slam area as well as a slightly higher pressure. The 
assumption is made that slamming will still occur forward of 
the area that is overestimated and that the pressure can be 
corrected by a relative angle correction. This will be 
discussed later. 
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SLAMMING PRESSURE AND HULL FORM COEFFICIENT (k) 
To estimate the slam pressure it iS necesSary to find a 
functional relationship for pressure vs. velocity. As mentioned 


earlier, it has been showi that equation (2) 


holds very well and can b?: used if the form coefficient, k, 
can be determined. 

Ochi and Motter [16] used mapping and regression analysis 
techniques to find the k values from experimental data. This 
method was used to develo» a computer program reported in 
Reference [16] Maecenas the coefficients. This program 
was used aS a Subroutine .n this work but it did not give 
satisfactory results for »road beam sections. The regression 
analysis had been perform:d over relatively narrow ship sec- 
tions and the confidence -ange for the coefficients was 0.073 
to 0.190 with 10% risk of error. In applying their program 
two problems developed: 

1. The k values for broad beam flat sections were 

well outside the 0.190 confidence figure; and 
2. The mapping blew up in many of the broad beam 
cases. 
Since many ships of comme-cial interest could not be covered 
by the program, another m:thod had to be found for these wide 


forms. 
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M.D. Ochi [17} in work which was an attempt to correlate 
the k values from two and three-dimensional drop tests as well 


as model tests proposed that the following holds: 


NO 


where B' half-breadth at 1/10 design draft 


ae 


half-section area at the 1/10 design draft 
This work, which was done prior to the regression analysis 
technique in Reference [16], was compared to the results ob- 
tained by the computer program. Better agreement inside the 
confidence range was obtained if a slightly higher constant 
was used. Therefore it is assumed that 


ono -— (6a) 


In the program developed in this thesis, equation (6a) is 
used when section values above 0.190 are obtained from the 
regression subroutine. This is done since equation (6a) yields 
larger k values and high pressures which is more conservative. 

Now that the form coefficients can be calculated, it is 
possible to determine the maximum pressure at each station 


along the ship if the velocity of impact can be found. 
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Mos. PROBABLE EXTREME PRESSURE 

The impact relative velocity 1S approximated by the 
amplitude of relative velocity which is a random variable 
following the Rayleigh probability distribution [14]. Using 
equation (2) and the Rayleigh distribution function, the 


slamming probability distribution function is given by 
f(p) = A exp [-A (p - p.)] (7) 
cone 
where p = slamming pressure 
ea threshold rressure = Kt? 
Ke= 1/ (Kk Ry) 
By applying order stitistics to the probability distri- 
bution (7), the distribution of extreme value of pressure in 
metinice number of impacts, Py? 1S given by 
n= il 


eee) tr (re) } Nesp 


£(p.) . 


ee ee) 


niexp [-A(p, - Po)] {l-exp[-A (p, 7 Po) ] 


Za 





ProlmedUdtion (8) 1t 1S possible to derive the most 
probable extreme pressure in n-observations, Po: Figure 3 
shows the probability distribution function given by equation 
(8). The most probable extreme pressure, Ps can be obtained 
by setting the first derivative with respect to P, equal to 


Zero. 


d £(p.) 


———— si = sé (2) 9 
‘ (9) 


Using equations (8) and (9), P is obtained 
+ Ro. ln(n) ) (10) 


This most probable pressure can be used to calculate slam 
forces and is one of two pressure options in the computer 


program. 


PeetkeMiE PRESSURE FOR DESIGN 

Although it may be reasonable to compare P, with the 
largest value observed in experiments, the probability that 
Pmemessure is higher than P, is 0.632. Since the probability 
Ge@serrbution function, equation (8) as plotted in Figure 3, 
1S never zero at high pressures, the extreme pressure for 
design can be determined only if a small probability of being 
exceeded (a) is allowable. The extreme pressure for design, 


Pa), can be found from 
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£(p_) d PL = el (11) 


eps aaa 


USing equations (8) and (ll) 


DB, (a) = k [r? = Re ee (e2)) 
The extreme pressure for design is controlled by the pre- 
asSigned small probabili’'y, a , which can be specified by the 


designer. In the calculations in the computer program Ochi's 


[14] value of x = 0.01 1S used, but can be changed if desired. 


SEA STATE INFORMATION 

In the calculations necessary to determine slam pressures 
the sea State inuSt be specified to determine 

1. the relative motions and velocities and 

2. the duration of the storm 
in order to calculate the number of slams from equation (4). 

In order to make the application of the program as simple 
aS possible, the only pa ameter which needs to be specified is 
the Significant wave hei: ht, H(1/3). This is made possible 
by the work done by Louk kis and Chryssostomidis [15] from 
which the r.m.s. relativ. motions and r.m.s. relative velocities 
may be determined for th: forward stations using only the 
Parameters of: 

meee = Block coef icient 


B 
2. L/B = Length to Beam ratio 
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So, = Beam to Draft ratio 

4. H(1/3)/L = S = Non-dimensional sea state 

5. F = Froude number = SHIP SPEED//gL 
iaese Darameters are sufficient to give good initial design 
values for relative motion and velocity. There is also a 
dependence on longitudinal center of flotation and radius of 
gyration but it 1S not particularly sensitive to these 
parameters [15]. These tables were developed from Series-60 
hull forms and, therefore, can only be used for ships which 
can be approximated by th2se forms. This allows calculation 
of slam forces for these ships without running seakeeping 
programs. The option is 1ilso available to specify the relative 
motions and velocities if Cesired. 

The storm duration éelso needs to be known. Following 
Ochi and Motter [14], wav2 records have been analyzed in the 
North Atlantic. A plot cf the envelope of storm duration vs. 
Significant wave height is shown in Figure 4. This envelope 
is used to determine the maximum duration to be expected for 
a ship operating at the significant wave height. This envelope 
is the default option in the program and is considered the ex- 
treme storm for design. If desired, a storm duration may be 


Specified in the input deta to the program 
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FIGURE 4- STORM DURATION vs. SIGNIFICANT WAVE HEIGHT 
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SLAM IMPACT REGION 

Using the time of storm from the previous section as well 
m@emene tabulated r.m.s. notions and velocities, equation (4) 
mooe, evaluated at each station considered to determine the 
number of slams during ar extreme storm. The number of slams 
1S greatest at the bow ard decreases rapidly as the midships 
crossection iS approachec. The sections where N(t) g es 
wesumea tO be the longitidinal extent of slamming. 

From model tests anc dropv tests [17] it has been found 
that the slam pressure is essentially dissipated at a depth 
of 1/10 of the design drz-ft, T'. Therefore the pressure needs 
to be integrated only over the forward portion of the hull 


where N(t) = 1 and below 1/10 design draft. 


SLAM FORCE AND IMPULSE 

The information iS row available to determine the most 
probable and extreme pressure for design PL and B (a) at the 
keel for each station as well as the slam impact area for a 
particular ship, speed, <nd Sea state. To integrate temporally 
and spatially, more infoimation is needed about the spatial 
variation of vressure in the girth direction as well as the 
time history of the slam. 


The local pressure iime history as used by Ochi [14] is 


shown in Figure 5. The curation of the slam pressure, th: is 
given by 
€, = 0.1 vE7520 = 4.4 x 10> VE ey) 
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FIGURE 5 - SLAM PRESSURE TIME HISTORY 
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which 1S an average for a 520 FT vessel whether it is in the 
flat bottom region or the bilge. Froude's scaling law is used 
for ships of different leigths. The slam pulse is assumed 
triangular with equal ris: and fall times. The duration of 
slam 1S one area where it is felt additional work needs to be 
done. In Appendix A an a.ternate equation for t, is considered 


aL 
and presented here as equition (14). 


Seo 100 7B! (seconds) (14) 


Cr 
i 


where B' = Half-beam a_ 1/10 design draft in feet 

Until further work can be done to experimentally determine 
which time duration is be;t, either can be used. The computer 
program uses equation (13 but it is felt by this author that 
from the analysis of drop test results beam should be used as 
the scaling parameter ins .ead of length. 

The work done by Och asSumes that there is no dependence 
on angle between the keel and wave surface, and that the worst 
slam case is when the slai pressure travels the full extent 
forward or aft at the slowest traveling velocity. It is felt 
that the wave slope has a large effect in reducing the slam 
forces predicted by Ochi's slam model. Therefore a new model 
of slam is proposed which accounts for the dependence on wave 
slope. 

Numerous experiments have shown that the slamming becomes 
maximum in regular waves ‘hen the ship encounters waves of 


length L on the order of 1.0 to 1.1 times the ship length, L, 


5) 





[12], [18]. It is assume? that a similar condition has been 
approximated in irregular seas leading to the worst slam 
condition. For the relat ve velocities to be maximum, the 
Ship will be pitching dow1 in the vicinity of zero pitch 
angle. For this to occur the slam will take place in the 
trough of the wave. Figu-e 6 depicts this model of slamming. 
The slam pressure will travel from the forward perpendicular 
aft and from the aft slam station forward meeting at what is 
called here the slam cent2r. Assuming the wave encountered 
has height equal to the significant wave height allows calcula- 
tion of the slope between the keel and the wave surface. This 
angular dependence will b2 used to modify the slam pressures 
from what would be expect2d for flat impact. 

Chuang [9] performed a series of drop tests with varying 
deadrise angles. It is ~roposed here that the reduction in 
pressure with deadrise arjle can be used for similar angles 
in the longitudinal direction. Figure 7 is taken from 
Reference [9] and shows @ Significant reduction in pressure 
with deadrise angle. The region between 0 and 3 degrees is 
where the compressibility of the air comes into effect. It 
is assumed this is not ag pronounced in the longitudinal 
direction due to the escéve of air in the transverse 
direction. The higher argles are almost linear when plotted 
on log, paper and the slope remains almost constant for all 
impact velocities. From the linear region of Figure 7 a 


correction equation for pressure was found as 


36 








WVIS JO IS3QOW- 9 3HNS!IS 


423U2> wes | 


| 
| 


ai 








_——y _ - 
¥ 

= ne 
$o= 


2 eo eee 


IN DEGREES 





Ww ee fan} 


| = 
hahed 
1 = = 
| ee a 
| iar i ot = 
ia i | - 
es : ine “ bad 
{ ; | | | “y 
el | eee 17 3 = 
. 7 | ie cs ae a 
jee be Sys Bee Sess feet © 
= sag << 
i ; ; 1 
ne , A | = 
; : 
I 


~» 


N 


: P 
Naa 
eS Se 
< ie 


' “ ~ a ee 
bane. io. ee ad ee gf Ont + Oe ir € 
So a“ a. ' ae ee 
Te so } —— een a b a 
ee ~, co 
oe — ne A pf ae seaierreraree ge 
Sn oN ee 
Se 
f sig ae peace Ais a ae IN culm ae mies) oy 
eae eee es Mes ik con Ae 
a CS (leet — c 
= Cas c> Cc se 
o => Fe =) 
wre 


AM PRESSURE vs, DEADRISE ANGLE 


UM 


- MAXI\I 


FIGURE 7 


38 











nN 


-0.12 |06| 


p,(o,8) = ey 10 (a) 
where Bp (o, 8) = maximum keel pressure with angular correction 
Pp (a) = pressure predicted by equation (12) 

6 = angle between the keel and wave surface 
in degrees 
If the wave is assumed approximately sinusoidal the wave 
slope can be calculated as 
6 = 935. H(1/3) X / L* (16) 


where X = distance away from the slam center in feet 

Using equations (15) and (16) the most probable or the 
extreme pressure for design may be modified to allow for other 
than zero impact angle. 

The peak pressures do not occur at the same time along 
the length of the ship. From full scale meaSurements on the 
Wolverine State and model tests [14], it was determined that 
the pressure pulse traveled from forward to aft or vice versa 
at speeds from 260 to 520 fps for a 520 ft ship. Ochi assumes 
that the slowest traveling velocity is the worst slam condition 
and that the results can be Froude scaled for different length 


ships to get 


Vers wa0ey 7520 = 11.4 7L ole) 


oS 





where a traveling velocity in feet per second 
This assumption was not felt valid since the traveling velocity 
should increase as the rate of entrance into the water in- 
creased--that is, as relative velocity increased so should the 
traveling velocity. It is agreed that if a slam could occur 
at the highest possible relative velocity and the slowest 
traveling velocity it would be the worst case, but it is felt 
that these two conditions cannot be met at the same time. 

The traveling velocity is calculated using the slam 
model in Figure 6. Using this model, the height above the 
water surface at the slam center, the most probable or extreme 
relative velocities of the slam center, and the length of the 


slam area are known or can be found. From this information 


the traveling velocity is given by 


Me 1i/2t, (17a) 


where 1 = slam length from FP to aft extent of slam 
ty = time slam takes to travel from aft extent 
of slam to slam center 
and 
t, = H/VEL (17b) 
where H = height of slam center above the wave surface 


at start of slam 


VEL = most probable or extreme relative velocity 
of the slam center 
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The last information required before integration can 
proceed is the spatial distribution of pressure along the 
girth. The most probable and the extreme pressure for design 
as calculated by egquations (10) and (12) and modified with 
equation (15) for wave slope are only valid at the keel. 
Following again Ochi's work [14], [18], the maximum slam 
pressure can be assumed constant along the flat bottom of the 


hull. From the end of the flat the pressure drops off approxi- 


mately linearly to zero at the 1/10 design draft, T'. Figure 
8a shows this for a typical section. The pressure component 
of interest is the vertical component on the hull. The values 


in Figure 8a must be multiplied by COS8 where 8 is defined as 
the angle between the normal and the vertical at the point in 
question. Figure 8b shows the new vertical component of the 
pressure distribution needed to find the slam force. 

All the information is now available to calculate the 

slam force and impulse. The procedure is as follows: 

1. Using eguations (4) and (5) along with the r.m.s. 
relative motions and velocities from tables [15], 
calculate the nunber of slams at each station. 

The last station aft with N(t) E 1 is the longi- 
tudinal extent of slam. 

2. Calculate the mo3t probable extreme pressure from 
eguation (10) or the extreme pressure for design 
from equation (12) at each station forward of the 


last slam station as found above. 
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FIGURE 8b - DISTRIBUTION OF VERTICAL PRESSURE 
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3. Assuming the certer of slam is midway between the 
FP and the last station where slamming occurs, 
calculate the slopes using equation (16) and 
the corrected pressures at each station using 
equation (15). 

4. Calculate the vertical pressure components at 
each station for approximately 10 waterlines between 
the baseline anc T' as in Figure 8. The spatial 
distribution of the vertical pressure is now known 
over the entire slam region. 

5. Assuming a triargular shape as in Figure 5, and 


a pulse duration of t, as predicted by either 


i 
equation (13) or (14), add temporally at each 
location. The start of slam pressure rise (t 5) 
at each spatial point is predicted by equation (17a) 
in the longitudinal direction and by the relative 
velocity in the girth direction. 

Due to the numerical integration procedure necesSary, a 
computer program was written to perform the above calculations 
and numerical integrations. The program and the details of 
the above procedure are discussed in the next section. 

By evaluating the pressures at different time increments, 
the force time history can be determined for each section. 


Integrating these force histories over time will give the 


slam impulse at each station. 
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V. COMPUTER PROGRAM AND DISCUSSION OF RESULTS 

A computer program was written to implement the procedure 
discussed in Section IV. A listing of the program is included 
as Appendix B. There are many input and output options which 
are discussed in Appendix C which is designed to be a user's 
manual. A short example program illustrating the data deck 
and some of the possible outputs is included as Appendix D. 

In this section the results obtained from the program will be 
discussed usSing the Mariner and Fotini-L as example ships. 

Throughout this thesis stations are numbered from 0 at 
the forward perpendicular to 20 at the after perpendicular. 
Since slamming occurs in the forward portion of the ship where 
the lines change rapidly, half-stations were introduced num- 
bered from 0 to 40 similarly. The slam forces calculated are 
considered as centered at the half-stations extending equal 
distances forward and aft of these half-stations. 

Two computer programs were written and used. The assump- 
tions and calculations of Reference [14] were duplicated and 
used to compare with the results in this work. The only 
deviations were the use of equation (6a) for k > 0.19 and the 
use of seakeeping information obtained from the Reference [15] 
tables. The results for the Mariner were identical to those 
obtained by Ochi [14]. The significant aspect of this is 
that good results were obtained without recourse to seakeeping 


Gonmpueer Programs. 
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OPERATING CONDITIONS 
The Mariner was used as one of the example ships and the 
computer program was run tinder various operating conditions. 


The four conditions selec:ed for discussion here are: 


1. MARINER FULLY LOADED (24 SKNOTS 
2. MARINER FULLY LOADED fS.> NOs 
3. MARINER LIGHT DRAFT 5 SK NOMS 
4, MARINER LIGHT DRAFT 7.4 KNOTS 


all at a Significant wave height of 25 feet. Table la lists 
eneeship Characteristics for the two loading conditions con- 
Sidered. This ship and conditions were chosen since they 


were used in Reference [li]. 


MARINER FULLY LOADED LIGHT DRAFT 
BLOCK COEF. (C,,) 0.613 0.583 
L/B 6.95 6.95 
B/T 2.82 3.80 
DRAFT FWD. 27 FEET 16.4 FEET 
DRAFT AFT. Dp 23.6 " 
H(1/3) 25" 25 : 
LBP 528" 528 4 
BEAM (B) 76 =O" 76 " 
DESIGN DRAFT Doe 7S DO Sees 


TABLE la - MARINER OPERA’ ING CONDITIONS AND DIMENSIONS 
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A second ship, the F)TINI-L, was analyzed and is one of 
the three ships reported ».n in Reference [4]. The principal 


characteristics and loadiig conditions are given in Table lb. 


Ee LENT =I, FULLY LOADED BALLAST 


BiGOCK COEF . CR 
L/B 
By T 

DRAFT FWD. 

DRAFT AFT. 

rae 3 

LEP 
BEAM 
DESIGN DRAFT 


SPEED (KNOTS) 


0.84 0.65 

Wile OD eS 

Cu SD 

Si 05 a ea 204 0b) B Ee 
5c, 0 Sl Cae 


20 to 40 FEET 


800 PEED 

106 " 
Sosume - 

Ooo ear oie 


Zee Ome Oe hE 
SUG EEED 
106 " 

Sie) 8) 


Se ae SC 


TABLE 1b - FOTINI-L OPERATING CONDITIONS AND DIMENSIONS 


R.M.S. RELATIVE MOTIONS AND VELOCITIES 


The first modification to the Ochi-Motter approach was 


to apply a broadness factor correction to the r.m.s. relative 


motion and velocity. 
narrow band processes, 


be obtained directly fron the tables. 


If the ship responses are assumed to be 
the r.m.s. values of the responses can 


These tabular results 


do not take into account the fact that seakeeping processes 


are not ideally narrow bend [15]. 


More realistic values can 


be obtained if a broadness factor correction is applied: 
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meen Gl = 6-2) (18) 





where m, = r.m.s. valve of response from tables 
ms = corrected r.m.s. value 
—€ = broadness factor 
and 
pa m2 
ao (Lo) 
o 4 
7 Gh 
where my = the k moment of the spectrum. 
Peevabue Of €¢ = 0.59 is used in the calculations [7], [15]. 


pes eCOrrection factor reduces the relative motions and 


velocities and therefore reduces the slam forces predicted. 


pemublo USING OCHI-MOTTEF METHOD 

The results for the four operating conditions using the 
Ochi-Motter method with the broadness factor correction are 
given in Tables 2 to 5 and Figures 9 to 12. These results are 
for each of the half-stations where slamming is predicted. 
Table 6 shows the total slam impulse and the amplitude and 
duration of an equivalent slam for the ship. Table 6 results 


are presented for ease of comparison of the two methods. 
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MARINER TOTAL SINE FORCE APPROX. TOTAL IMPULSE 


CONDITION AMPLITUDE DURATION (KIP-SEC) 
(KIPS) (SEC) 
1 985 0.43 269 
2 4490 0.62 1766 
3 4406 0.68 1901 
4 1220 0.48 372 


TABLE 6 - TOTAL SLAM FORCE AND IMPULSE (OCHI METHOD) 


PEW METHOD RESULTS 

The changes made to the Ochi method discussed in the 
previous sections and incorporated into the second computer 
program are: 

1. an angular correction applied as per equation (15), 

2. a different model of slamming as shown in Figure 6, 

and 

3. a traveling velocity calculated by equation (17a). 
The results for the same four operating conditions are given 
in Tables 7 to 10 and Figures 13 to 16. Table 11 gives the 
total values of slam force and impulse corresponding to 


Table 6 above. 
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MARINER TOTAL SINE FORCE APPROX. TOTAL IMPULSE 


CONDITION AMPLITUDE DURATION (\KEPSeHe) 
ChE? sS) (SEC) 
el G2 Ora 3 rs 7 
2 aoa? Ora2 1 I78 
5 O'S Gea 2 oo) 
4 2084 OF 7 22 


TABLE ll - TOTAL SLA‘ FORCE AND IMPULSE (NEW MODEL) 


Seen SON OF RESULTS BETVEEN THE TWO METHODS 

In general the new model predicts lower slam impulses, 
shorter slam durations ani higher slam forces. [In addition, 
the time histories are moce triangular in the new model. 

The lower slam impulses are a result of the angular 
correction term only. If the correction were not used, both 
methods would yield the same impulse values. 

The shorter total duration is a result of the slam model 
and the calculated traveling velocity. The model reduces the 
duration by approximately 1/2 since the pressure wave travels 
forward and aft toward the slam center. An additional decrease 
in slam duration occurs because the calculated traveling 
velocities are, in general, higher than those calculated by 
Ochi. For these four operating conditions the traveling 
velocities range from 731 fps for case 3 to ]016 fps for 


case l. 
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The total force on the ship is the spatial and temporal 

integration of the local pressures. In the new method the 
pressures are always less than or equal to the Ochi method 
but the new method yields larger forces. This is due to the 
mitegration Over time. Referring to the plots of both methods, 
Figures 10 and 14, the reason for this can be seen. The half- 
Station forces are more Eunched together in the new method. 
The total force on the ship at any time is the sum of the 
menpces On C€ach half-station at the same time. Due to the 
greater overlap in the new method, more half-staticns are 
summed at the same time, resulting in higher forces. 

The importance of the various parameters of storm duration, 
threshold velocity, local slam duration, significant wave 
height, and most probable vs. extreme pressure for design 
were looked at to determine the sensitivity of the results 


@omchneir Variation. 


SENSITIVITY TO STORM DURATION 

To determine how important accurate prediction of the 
maximum storm duration was, the MARINER case 2 was run for 
fiererent Values of storm duration with all other input 
variables held constant. AS Seen in Figure 17 the total slam 
impulse is not sensitive to storm duration. A 50% increase 
in storm duration from 3C to 45 hours results only ina 4% 


increase in slam impulse. 
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SerolLlivViTyY TO THRESHOLD VELOCITY 

Since the threshold velocity at which slamming occurs 
has been defined differently by various authors, the sensi- 
tivity of slam impulse to changes in threshold velocity was 
determined. Values of threshold velocity between 6 and 20 
ft/sec were used for MARINER case 2 and the total slam impulse 
plotted against threshold velocity in Figure 18. From this 
figure it can be seen that the use of 12 ft/sec is reasonable 


and, at worst, conservative. 


SENSITIVITY TO LOCAL SLAM DURATION 

The total slam impulse on the ship is the sum of all the 
local impulses. Those parameters which affect local impulse 
also influence the total impulse the same amount. From 
Figure 5 the local slam impulse is given by: 


eee? 2 Pp (a) ¢t (20) 


iu 1 


where I, = local eee ene = sec ania 

From (20) it can be seen that a large variation in ty 
results in a large change in I, and subSequently in the total 
impulse. Therefore tota] slam impulse is very sensitive to the 
local slam duration. Sirce the determination of ty is in 


question, this is considered the area where further work is 


necessary (Appendix A). 
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AFFECT OF SIGNIFICANT WAVE HEIGHT 


Another ship, 
computer runs were made uSing 
ficant wave height from 12 to 
The Ochi method 


the results. 


times the new model results. 


the Fotini- 


ereases with 


seaway. The 


Significant wave 


Fate Of increase 


L, waS also analyzed. Several 
both programs varying the Signi- 
40 feet. Figure 19 is a plot of 
gives results approximately 2 1/2 
As expected, the impulse in- 
height and with velocity in the 


decreases at higher Significant 


wave heights and is expected since the ship responses behave 


Similarly at higher sea states. 


MOST PROBABLE vs. EXTREME PRESSURE FOR DESIGN 
The four Mariner cases were run using the new method for 
the most probable extreme pressure and the extreme pressure 


for design cases. Tables 11 and 12 summarize the results. 


MARINER TOTAL SINE FORCE APPROX. TOTAL IMPULSE 
CONDITION AMPLITUDE DURATION (fae She) 
hones) (SEG 
1 USS OtZ.0 it 
Z =248 2,5 510 
3 OO. OGe2 > Bo 
4 1147 OFS 137 


TABLE 12 - TOTZL SLAM FOHCE AND IMPULSE (MOST PROBABLE 
EXTREME PRESSURE USED) 
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The duration of the slam is greater in the most probable case 
Since a lower relative velocity predicted results in a slower 
eraveling velocity (n.b.,Table 11). The impulses are approxi- 
mately 1/2 those obtained from the extreme pressure for design 
calculations. Since the most probable extreme value is indi- 
Ccative of the largest value observed in a trial, these values 


should be used if any comparison is done with full scale data. 


SHAPE OF THE SLAM FORCING FUNCTION 

Knowledge of the shape of the slam forcing function is 
necessary to input the slam force into the S.H.V.R.S. program. 
The S.H.V.R.S. program will accept three input shapes for the 
poromng function: 

1. triangular 

2. half-sine wave 

3. damped half-sine wave. 
The forcing function can be applied at each of the half-stations 
with time delays corresponding to the traveling velocity or a 
Single slam force can be applied at one location. Since much 
of the previous work [2], [3], [4] used a single sine pulse, 
this program develops a half-sine approximation to compare with 
these results. 

The shape of the force at each half-station is different 
in the Ochi method and this method. In the Ochi method the 
assumed shape of the local slam (Fiqure 5) is not critical. 


Several pulses of different shapes were tried in the Ochi 


WZ 








program but the predicted half-station forcing functions were 
about the same. This is because the traveling velocity is 
Slow enough that the slam starts at different times along the 
length of the half-station. In the new method the traveling 
velocity is, in general, high enough that the slam starts at 
the same time along the entire half-section. This results in 
a force shape the same as the assumed local slam pulse. The 
shape of the slam force at each half-station (Figures 9 to 12) 
in the Ochi method can be approximated as half-sine functions. 
The program calculates the amplitude and duration of these 
approximations. The shape of the force at each half-station 
(Figures 13 to 16) in the new method can be assumed triangular, 


of height equal to the maximum force, and of duration 


t) = 1.1 t, = 4.8 x 100° VL (21) 


where t. was determined by neglecting the "tails" of the 
force function and maintaining the same impulse. 

The total slam force shape can be obtained by adding the 
forces from each half-station at the same time. Figure 20 is 
a plot of the total slam force for the MARINER case 2, time of 
storm 30 hours, calculated by the new program. Figure 21 is 
the plot of the forcing function at each half-station. If the 
points which are less than 1% of the maximum value are neglected 
and a linear regression analysis done on the remaining points 


in Figure 20, the triangle formed is almost equilateral. There- 


fore this shape should be used for the total slam forcing »fEuncecion. 


12 





12 


10 


(OL X SdIW) 3DYOS WVIS WW1LOL 


0.2 


0.1 


TIME (SEC) 


F 


74 





SYNOY OLE =WHOLSL ~- TSCOW MAN - YSNIYVW-S - te 3YNSIS 


(335) SAIL CSNIVAN 
$*Oo ¥ “Oo E“O 2°0 
Mg ol 


Ny 


Wy 

















SEERERER BRO 


(SdI¥) SOu04 





MIDSHIP SLAM BENDING MOMENT AND STRESS 

The FOTINI-L was considered at different significant 
wave heights, velocities, and loading conditions. The results 
for the total slam impulse in the ballast condition are shown 
in Figure 19. As the significant wave height and/or ship 
speed is increased, the slam center moves aft. Donovan [4] 
determined midship bending moment and stress for the FOTINI-L 
using a 100 ton-sec half-sine slam impulse, applied 36 feet 
aft of the forward perpendicular. This slam location is close 
to half-station 2 of this model. From the computer runs on the 


FOTINI-L, two operating conditions resulted in slam centers at 


half-station 2. These are given in Table 13. 
CONDITION LOADING SPEED H (173) TOTAL DURATION 
CONDITION (KNOTS) (FEET) IMPULSE OF IMPULSE 
(KIP-SEC) (SEC) 
-POTINI-L FULL 14.38 32 ie OES 
2-FOTINI-L BALLAST oro 20 960 Oss 


TABLE 13 - FOTINI-L ~ SLAMS PREDICTED AT HALF-STATION 2 


From the graphs of response in Donovan's work, these 
slam impulses give the moments and stresses in Table 14 for 
various stiffnesses. This is determined assuming a linear 
relationship between bending moment and impulse which holds 
if the slam durations are the same--that is, holding all else 
constant, doubling the slam impulse doubles the bending 


moment response. 
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DELP PNESS* 


CONDITION 60% 80% 100% 120% 140% 
BENDING BE OriNnit—1, 343 364 375 391 404 
MOMENT 
fe r-TONS ) 2-FOTINI-L 246 268 283 296 SINS) 
STRESS L=Pe rir 1, AO Dee eo EG. > 14.5 
(KPST) 


Zoe TN eal, Zee nee Ge aD Loree 


TABLE 14 - FOTINI-L ~- BENDING MOMENT AND STRESS AMIDSHIPS DUE 
TO SLAMS CENTERED AT HALF-STATION 2 


*The stiffnesses are % of as built values. 
The values in Table 14 for the fully loaded case are 


approximately 5.4 times what was predicted in Reference [4] 


and 4.3 times the predicted results in the ballast condition. 
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VI. RECOMMENDATIONS FOR FUTURE WORK 

AS seen by the values of stress obtained in Table 14, 
the slam forces and impulses result in large midship stresses 
on the FOTINI-L. The stresses could be evaluated Only so 
slams at half-station 2 since dynamic response results were 
available only at this station. As the slam center moves aft, 
the total slam impulse and force increases. Higher slam 
Stresses could result from slams at different locations. 
Therefore additional slam locations should be considered in 
the S.H.V.R.S. program in any study of the slam response. 

The shape of the total slam force is more triangular 
than sine (n.b., Figure 20). The S.H.V.R.S. program should 
be run for several triangular pulses and the results compared 
to the half-sine results available. If they do not come close, 
a triangular pulse should be used. 

The validity of using a single total force approximation 
vice distributed forces along the bottom needs to be verified. 
The results of this work give the distributed forces as well 
as a half-sine approximation for the total slam force. Both 
Should be run on the S.H.V.R.S. program to determine if the 
two agree. If they do not, either the distributed forces 
should be used or a new approximation developed which does 
agree. 

Although the slam impulses are less in this method than 
in the Ochi method, they still result in large slam stresses. 


There is still reason to believe, because of these large 
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values, that these impulse values are also overpredicted. As 
@iscussed earlier, and in Appendix A, equation (13) for the 
duration of the local slam pulse is in question. There pre- 
sently exists a large amount of literature on drop tests but 
few model and full scale results are available to determine 
durations. Correlation with full scale data needs to be 
accomplished to determine if equation (13) is valid or some- 
thing like equation (14) would be better. 

A series of model tests should be run to verify that the 
model of slamming used here is valid. This would include 
verifying or modifying the traveling velocity calculations 
and the correction equation for the angle between the keel 


and the wave surface (equation (15)). 


eo 


wil. CONCLUSION 

The work presented here develops a slam force prediction 
method which is sensitive to seakeeping parameters based on 
a model of slamming different from other methods developed 
to date. The results obtained show slam impulses and time 
durations less than those in the Ochi-Motter method [14], but 
impulses greater than the 100 ton-sec value used by Kline [3]. 
It 1s felt that values closer to those predicted here should 
be used in the overall A.B.S. study on hull stiffness. Asa 
minimum, the results of this work have shown that the 100 ton- 
sec value used in the work to date is in question and should 


not be applied to all ships under every operating condition. 
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APPENDIX A - LOCAL SLAM PRESSURE, DURATION OF PULSE 
Ochi predicts that the slam pulse is triangular and the 
duration is 0.1 seconds for a 520 foot ship. From this he 


Froude sacles and concludes that equation (13) holds - i.e., 


t, = 4.4 x ont 

From full scale measurements on the USCGC UNIMAC [19], 
L = 300 feet, values of slam duration locally were of the 
order of 0.020 to 0.025 seconds. USing equation (13), this 
length ship results in a predicted duration of t, = 0.076 
seconds. This discrepancy was too great to ignore without 
further consideration. 

From drop tests done by Chuang [8], [9], and [20], on 


flat bottom impact with zero deadrise angle, the following 


equation holds for deadrise angles less than 3° 


t) = 4 ey ee (22a) 
where 

LS wee wicth Of infinitely long flat plate 

a Speed of sound in air = 1125 ft/sec 


For the UNIMAC and the durations of slams measured this 
would predict the half-width at frame 23 to be of the order 


of 6 feet. Lines were not available to determine how close 


8 3 





this agrees to the half~-breadth of the I7/VOsdratt dt ermirs 
section, but it is a reasonable value for a ship with beam 
ee tt. 

In other drop test data [21] equation (22a) held quite 
well for many different sized test specimens if the beam was 
used as the scaling parameter. Without more information from 
full scale tests, equation (22a) cannot be verified and there- 
fore will not be used. It is felt that the local half-breadth 
Should be the scaling parameter, not the length as used by 
OCT . 

It 1S proposed here that 


3 


= 3.6 x 10 ~ B! (14) 


a 
which comes from equation (22a) would be a better estimate 
of slam duration with possible modification of the constant 
from full scale data. 

Equation (22a) holds for small deadrise angles less than 
3° where the compressibility of air influences the duration of 
pressure pulse. For larger angles the impact pressure time 
history becomes more steep fronted with about half the duration 
predicted by equations (22a) and (14). Therefore above 3° 


equation (22b) holds fairly well. 


_ 22b 
ty 2 eG ( ) 
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Although equations (22a) and (22b) are for deadrise angles, a 
Similar result might be expected for angles between the keel 
and the wave surface. Therefore equation (22a) is probably 
valid in the flat region near the slam center while equation 
(22b) 1s probably valid outside this region. 

For the MARINER equation (14) yields values of t, = 0.022 


i 


at half-station 4 to ty = 0.115 at half-station 16. Therefore 


it results in much lower impulses than would result using 


equation (13). For the FOTINI-L larger values were found 
uSing equation (14) than from equation (13). At half-stations 
4 and 10, ty = 0.130 and 0.189 respectively using (14). This 


would result in larger inpulses than one would get if (13) 
were used. If equation (22b) were applied in the appropriate 
regions, the above resulting impulses could be reduced. 

Use of equations (22a), (22b), and (14) are not justified 
at this time, but it is felt that full scale data should be 
analyzed to determine if equations of this form are valid. If 
they are, further reduction of calculated impulses by the 


program is possible in many ship cases. 
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APPENDIX B - COMPUTER PROGRAM LISTING 

This appendix contains a listing of the FORTRAN computer 
program. It was written for an Interdata Model 70 computer 
and requires approximately 60 K of memory. Two subroutines 
are called which are not listed here but are readily available 


to most users as part of the IBM Scientific Subroutine Package. 
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PROGRAM SLAMCW WALKER JANUARY 1976 
PROGRAM MAIN -== THIS PROGRAM CALCULATES THE StaM FORCES AND 
IMPULSES GENERATED ON &A SHIP IN AN IRREGULAR SEA 
DIMENSION XIN(S) sYIN(S) sXOUT(11) eYOUT (11) oxuCl7oll)oVYY(1lL1) eT(20)% 
c VK (17) eYOC12) eQ®M617) oPV6(17) eEPRES (17) ePRNAPR(17), 
Cc PRES(17011) eFORCE (17220) sNVEL (17) e TIME (17) 
DIMENSION A(17) +DUR(17) eSLAMS(17) eTM(20) eSHAPF (17) 6 
CPRSLAM(17) eFMAX(17) eAIMPLS (17) oPVEL (17) eXGC17011) oe XSCL (4) 
DATA ESL eESPeNINoNOUT 0 P1799 6099 0 0S0 11930 141592654/ 
DATA NeM/8.,S/ 
DATA SHAPE /0 89375 20.791] 6667 40.69375 00.690.5106167 60 0425 0% 34759 
CO 6266666 0 0.1937500.125 506 1604166620. 000-02056254-0.1083332°0.15625 
C2=02.20-0.239583/ 
} REAN (Nel000) 101 +102-104e108 
1000 FORMAT (71S) 
IF (102.LT.1) 60 TO 999 
Cc READ SHIP NAME AND LOADING CONDITION 
READ (No1001) (TCI) ¢ 121-20) 
‘ 3002 FORMAT (20A4) 
c READ SHIP DATA CONSTANT IN PROBLEM 
REAN (Ne1002) RR eHHe ALBPeALWL e STASP 
1002 FORMAT (SF19,5) 
WRITE (Me1003) : 
1003 FORMAT (°1%+35xXe*THIS PROGRAM CALCULATES THE SLAM FORCE DISTRIBUTI 
CON VS TIME %./s38Xe*0N THE LOWER 1/10 OF THE SHIP FROM STATION 8 FO 
CRWARD.%e////) 
c WRITE SHI& NAME AND DATA 
WRITE (Me1lNNG) (TCI) eo 121920) eBBeHHe ALBPoALWLe STASP 
100% FORMAT (20Xe2084e//e30X—e"HALF BEAM = %oF1N.S5e10Xe*DESIGN ORAFT = * 
CoFlLOeSe// oe 30Xe (LRP = %oF1062016Xe*LENGTH LOAD WATFR LINE = *2F10.2 
Co//e30Xe DISTANCE BETWEEN STATIONS = *%oF Ge2e40Xe//e35Xe* (ALL DUESTA 
CNCES IN FEET) *s////) 
Cc CALCULATE THE K VALUES USING SUSROUTINE SPOON 
i H=HH/10, 
S=HH/100- 
XOUT(1)=0.0 
00 10 I[=2611 
XOUT (TIT) =XOuUT(Te1) ¢ S 
10 YY(I) =xXOUT(TI) 
re I=z0 
Cc READ OFFSETS AT ANY S WATERLINES 
READ (Nel005) (XIN(J) eJel eS) 
100S FORMAT (5F10.5) 
20 READ (Ne1l005) CYINGY) »Jel eS) 
Iz=Ie} 
IF (YIN(S).LE9-5S) GO TO 30 


‘eke ka) 


CALCULATE OFFSETS IN 10 INCREMENTS FROM Bele TO 1/10 DESIGN ORAFT 
USING SUBROUTINE UGLYDK 


AAA 


CALL UGLYDK (NINoNOUToXINe VINe XOUT o YOUT sESLeESR) 
00 21 J=1el1) 
JJ=ZIL2<J 
YO( J) =YOUT (JJ) 

2] xXx LT eJ) =YOUT (J) 
B=YOUT (II) 
TaOP=I01 

Cc CALCULATE THE FORM COEF. OF THIS STATION 

CALL SPOO (YOoVKK eB eHeSe3B eHHe LA0P oI) 
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VK (IT) =VKK 

IF (104.LT.1) GO TO 20 

O=FLOAT(I)/2. 

WRITE (M1006) Oo (VYING J) oXIN( J) » YOUT (J) eXOUT (J) o S21 05) 9 (YOUT (J) 6 
XOUT (J) eo J=601)) 


Cc 
1006 FORMAT (91% 649% 6511999) pS eSOKe PHM ep GAKe PHP o/eSOKe POH eo LIXe PSTATION' 


30 
4003 


1011 


100 


1007 


1920 
1008 
4000 


500 


501 


110 


CFSele® OFFSETS (FEET) & > DK et Ot of eo SOK e POP e4AKe MFO yf ye SOK GSO (HHH) oS e 
C5OX so PFO GR MOH gfe SOK PA Pe GKe PF INPUT 8 oe IX eo 990, 7X *CALCULATEN% s 7X, 2%? 
CoF SOK o OFF PINK y MMH ye PUK oe HAM eS gp SOK SD (OFF) oS oe SOXe POH e lL 1K oH ge 1 DXQ PO? 
Col lUX se PFO oe LAK PAH gS ge SNK oh SK oO KH FSX eo PHOS oe MVM SH go PHOS aI XO SK Gg? 
CEH SX eo PV 6X9 MOH Ag HOOK HO( 9H) pS eSOKe PF Mo LIX oP HO eo LIX et HO ot Xe 8 O%y) 
COK a PFO gf oF (SOK e PH eo FR Py IK o POOF Bee Xo POO ge FA Pe 3K oO MOF BH PoGXe te? 
Co) ob (SOKe MH QoL 1X oP FO eo] LX 9 OHH FHA Ao Ko POH eA Lo GXs POH g/) oSOXCS0(9* 
C*)o/e*0*) 

GO TO 20 

NLINES=I-1 

READ (Ne4003) IPLT + (XSCLIK) »K=1¢4) 

FORMAT (1564F10.2) 

CALCULATE AN AVERAGE FORM COEF. FOR FORWARD PORTION OF SHIP 
CALL KEFFEC (VK-.AVGK) 

WRITE (MelO11) AVGK 

FORMAT (91% 639K 950 (998) pS eUOKe PHO ge GAX oP OY g/g GXy POC e 
C 6X0 * AVERAGE HULL SHAPE FACTOR = %eF Bee 6Xo POP oSy 
C4OK eo FF op GAK eg PHO eS eKOK 5019 F9) OSS) 

READ IN INTEGESS FOR CONTROL OF OUTPET OPTIONS 

READ (Nel00N) 1O09eT03e10Se106e107-1910 oe, 
IF (I09.LT.1) 60 TO 1 
READ (Ne1l007) (T( 1) eI=1-20) 

FORMAT (2944) 

NSTA=NLINES 

WRITE (Me1020) (TCT) eFT=1920) 

FORMAT (*1% 9 1Xe2MAGe//) 

READ (N-e1008) FROUDE + SPEED> SEAsHF eHAeCBe VAFT»VTHR 

FORMAT (8F19.5) 

READ (Ne%000) TSTORM 

FORMAT (F19.5) 

HSIG=SEA*®ALBP 

IF (TSTORM.GT.0.01) GO TO S01 

IF (HSIG.LE.35.) GO TO 5$00 

TSTORN=50.0-1.0°HSIG 

GO TO S01 

TSTORM=78.0-]1 -8*HSIG 

IF (HSIG-eGE.20.) GO To 501 

TSTORM=45.-0.15*HSIG 

CONT INUF 'e 
IF (FROUDE.GT.0.00001) GO TO 110 . ‘ 
FROUDE = SPFEN®1].6889/SQRT(32.174°AL WL) 

SPEED = FROUDE*SORT(32.174*ALWL) /1-6889 

FSCALE=SORT(ALAP/S20.) : : 

THETA = (HA-HF) /ALBP 

VAFTS=VAF T®F SCALE 

VTHRS=VTHR*F SCALE 

ALAB=ALBP/ (?.*88) 

ABT=4.°RB/ (KA*HF ) 

IF (105.€9.1) GO TO 103 

READ (Nel908) (RM(T) eRVI IT) o 1=228e6) 

CALCULATE THE RELATIVE MOTIONS AND VELOCITIES AT ALL HALF-STATIONS 
CALL RELMV (RM,RV) 

GO TO 104 
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103 NR=0 
lee IF (NR.GE-NSTA) GO TO 104 
READ (N»1008) RM(NR*1) eRV(NR*1]1) sRM(NR*2) oRV(NRO2) sRM(NRO3) 9 
CRV (NR 3) sRM (NRG) sRV(NRO4S) 
NRENROG 
GO TO 1l22 
104 CONT INUE 
Cc CALCULATE SLAM PRESSURES AT THE KEEL» NUMBER OF SLAMSs EXTENT OF 
c SLAMs VELOCITIES OF SLAMs ANDO PROBABILITY OF SLAM 
CALL SLAMPR (RMsRVoHF eHA SALWL »NSTA cE PRES + VK 9 SLAMS « PROBPReNVEL e VTHR 
CSeTSTORMsHSIGePRSLAMsPVFEL oSTASP) 
NCENTR = NSTAS? @] 
IF (VAFT.GT.9.01) GO TO 105 
IF (NSTA.FO.0) GO TO 195 
HI=(L.-SIN( (PI/2.)-0.0714*FLOAT (NSTA))) *HSIG/2,. 
IF (I106.EQ.0) GO TO 106 
TT=HI/PVEL (NCENTR) 
GO TO 107 
106 TT=H1/0VElL. (NCENTR) 
107 VAF TS=FLOAT  (NSTA) ®STASP/(4,*TT) 
VAFT=VAFTS/FSCALE 
10S IF (NSTAEQN.0.AND.IO3.EN-1) READ (Nel089) YoU 
1089 FORMAT (13) 
WRITE (M9114) CBeFROUDE sSPEED sALB ec ABT eSEAsHSI Ge TSTORMs THETAsHF op HA 
CeNLINES »NSTA 
1014 FORMAT (°0%.19Xe"THE SKIP OPERATING CONDITIONS ARF 2*o//> 
C20Xe*BLOCK COFF. = %eF 7, 40/5 
C20Xe*FROUOFE NO] = %oF7,.4925Ke "SHIP VELOCITY = oF G.20% KNOTS# e/g 
C20Xe*LENGTH TO BEAM PATIO = "oF 6.39/09 
C20Xe*BEAM TO ORAFT RATIO = %eFHe30/ 9 
C20Xe *NON-OIMENSIONAL SEA STATF = "oF 7TeGelOXe*SIGNIFICANT WAVE HEIG 
CHT = *eFS.00" FEET eS 
C2OXe* TIME OF STORM = #,F5e00% HOURS?# o/s 
C20Xe* TRIM ANGLE (RAOIANS) = %oF TeSe/o 
C20Xe* DRAFT FORWARD = % oF 5e2e24Xe* DRAFT AFT = oF S.20/¢8 
C20Xe*NUMBER OF SHIPS LINES READ = *,13913Xe*NUMBER OF HALF STATION 
CS CONSIDERED = *.1307) 
WRITE (MelN15) VAFTsVAFTSeVTHR ag VTHRS 
1011S FORMAT (20xX+*TRAVELING VELOCITY = *.F9.1le* FT/SEC FOR S20 FT SHIP 
C*e//o 
C2OXse* TRAVELING VELOCITY = *oF 9.le* FISSEC FOR THIS SHIP*//>5 
C2OXe* THRESHOLD VELOCITY = "ef 6.20" FT/SEC FOR S20 FT SHIPY ose 
C20Xe*THRESHOLO VELOCITY = "oF 6.le* FT/SEC FOR THIS SHIP* o//s 
C20Xe*THE PROBAAILITY OF EXCEEDING THESE SLAM FORCFS (ALPHA) = 0.01 
Cte//) 
IF (NSTA.EQ.0) WRITE (M,3030) 
IF (NSTA.FOQ.0) GO TO 19 
3030 FORMAT (40xX%e*NO SLAM CALCULATIONS PERFORMED SINCE THE % 9/7 e40Xo 
C*NUMBER OF SLAMS IS LFSS THAN 1 FOR ALL STATIONS'»//) 
IF (106.EQ.0) GO TO 150 
c CALCULATE SPATIAL OISTRIBUTION OF PRESSURE ALONG GIRTH 
CALL PRESS (PROBPR INSTA cPRES oO XX eo YY¥ 9S 0XG) 
60 TO 63 
150 CALL PRESS (EPRESeNSTS oPRES oe XXe VV 9S eXG) 
63 CONTINUE 
Cc INTEGRATE PRESSURE TO GFT FORCES FOR EACH HALF=STATION 


CALL FORCES (FORCE »NSTAc TIME eDELTAT o PRES eo XGeALBP oNVEL s VAFTS+STASP > 
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Cc YY ePVELe 106) 

CALCULATE HALF-SINE APPROXIMATION WITH EQUAL IMPULSE AND QURATION 
AT EACH HALF -STATION 

CALL APPROX ( NSTAsFORCE oF MAX sAIMPLSe AeDURsDELTAT) 

IF (1010-EQ.2) GO TO 3020 

IF (TOLO.FN.1) WRITE(M.2012) 

BOL] FORMAT (CPL L2K (EME) oS ot KXo2( POO IK) oe POU LGKg COOL TK eo MOH PG (LIX, HOt) 

CoFolXo2(*F9,9Kdo9%#,¢ REL MOTION *#® EXT PRESSURE 611% 


C#® MAXIMUM ®¢¢11X09® AMPLITUNE ®#,11X—"% TIME i A) 
CiXse* HALF ® FORM %934( 90), NUMBER # SLAM ° SLA 
CM * oF SINE ® SLAM . SLAM P%o/e 


CIXe**® STATION ® COEF * REL VELOCITY * MOST PROS PRESS ® OF SLA 
CMS ® FORCE * IMPULSE ®* APPROX ® QURATION ® STARTS @t, 
CLolXetFeoIx, 98 («) % (NON-0IM) PO e6XeF (PST) Pe6Xeot Feel iXes 
cee (KIPS) * (KIP-SEC) *® (KIPS) #e,e(e (SEC) PO) oS 
CiXsl26(*4e)) 

IF (TO10.E€Q.1) WRITE (M-3013) (JeVK (J) eRM( J) oEPRES (J) oe SLAMS(J) » 
CFMAX (J) eAIMPLS (J) 0 A(S) eMUR( J) e TIME (J) eRV (J) ePROAPR( J) »-J21eNSTA) 
JOL3 FORMAT (17(1X 09 F% eo 3BKoTBIe3Xe tM to AX oF Gu do® BH ge uXeF Be teGXet Bl yOXoF Hels 

CSXe 
COFHO eA(LXoFD.Ool Xo tl) oeAKeFRelolXe Pt eeXeF BelolXet Otel (3XoF be te2Xe 
CHPE) oS oT Ko Fe A(GK oF) pK oF He BeGXotB* oGX oF Gel oSXo tM QeG(LiXetHe)y» 
CJolXel26(9!) 7)) 
IF (1010-E9.0) 69 TO 3015 
IF (106.E9.90) WRITE (4.3011) 
IF (3106.EQ.1) WRITE (M.3014) 
3011 FORMAT (30X-°°#% EXTREME PRESSURE USED IN FORCE CALCULATIONS #e0e, 
Cvs) 
3014 FORMAT (30x-**8®% MOST PROBABLE PRESSURE USED IN FORCE CALCULATIONS 
Cc eee, s/) 
GO TO 3020 
3015 900 3000 KM=)] +NSTA 
TM(L)=TIME (KM) 
00 3006 KN=2920 
3006 TM(KN)= TM(KN=]) * DELTAT 
WRITE (Ms3001) «<M 
If (107.E0.0) GO TO 3010 
IF (106.FQ.9) WRITE (4.3002) 
IF (106.FQ.1) WRITE (M.3004) 
WRITE (Me3005) (TM(S) oFORCE (KMoJ) o J21 920) * 
WRITE (403007) 
2010 CONTINUF 
WRITE (Me3008) EPRES{( KM) sPROBPRI(KM) » VK(KM) 
WRITE (Me30N9) RM(KM) eRVIKM) p SLAMS( KY) ePRSLAM(KM) -OVEL (KM) ¢ 
Cc PVFEL (KM) oF MAX (KM) eAIMPLS (KM) 0 A(KM) eOUR(KM) » TIME (KM) 
3009 FORMAT (30x-.°NON-DIMENSIONAL RELATIVE MOTION = % of Te4e/o 
C30Xe*NON-DIMENSIONAL RELATIVE VELOCITY = *%oF Te4esy 
C30OXe*NUMBER OF SL4MS NURING EXTREME STORM = %oF B.oZe 
C3OXe*PROBARILITY OF SLAM IMPACT = %eF 8.5076 
C3IOX+*EXTREMF RELATIVE VFLOCITY = "oF Beleo® FISSEC* w/e 
C3IOX~e*MOST PROBABLE RELATIVE VFLECITY = oF Bale* FTSSEC* oS /o 
C3OXe "MAXIMUM VALUE OF SLAM FORCE = %oF 10ele* KIPSt o/s 
C3IOX+* IMPULSE VALUE = %oF 10.00% KIP*SEC# o/s 
C3IOXs* AMPLITUDE OF SINF APPROXIMATION = *oF 10ele* KIPS*e/o 
C30Xse*DURATION OF SLAM IMPACT FORCE = %oF 8.50% SEC ose 
C30OXe* TIME SLAM STARTS = %oF 8.50" SEC*o//) 
IF (106.E0.0) wRITE (4.30112) 
IF (106.F0.1) WRITE (4,3014) 
3000 CONTINUE 
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1 


1 


1 


c 


3020 CONTINUE 

3001 FORMAT (91% 629K e40( 8 OF) So 30K PHF IBXe OHH e/a BOX OO eTOXe HALF STAT 
CION Feldol [Xe FtoSog 3K e (Hho JAK oH FH So BOK HO 99%) ) 

3002 FORMAT (30K et Ft ol 7X et Ft yg IKe* EXTREME FORCE % BX oP Ft oS eo BNKX oH Fg BKO IT 
CIME (SEC) Ho 3Xo 88 eSXe MFOR DESIGNE SX oP St oS eg IOX eo MFQ oT IX 4 POH TK Ot (K 
CIPS) Pe TK oF FO gfe IOX e406 94) ofa BOX e POU eT TX et HE SHK, OOP) 

3004 FORMAT (30K ef Oto TTX eo Ft yg IXe*MNST PROBABLE %, 3X ot Ft oS eo ZOX ot HHA IKe IT 
CIME (SEC) © eo 3X0 Oho SXe! FORCE CoS K eM OMe Se TNX eo POO eo LT Xe Hho 7Xe ht (K 
CIPS) Ho TX el FH oS go DOK eG (988) oeSe BOX at Ft ol TKXo FH ge AOK, 9 HP) 

3005 FORMAT (39K oOo SK oF 7. seo tat ep AX oF Bel eGXe 0% ye) 

3007 FORMAT (30% .49(°%8!) //) 

3008 FORMAT (30XetEXTREM“E POESSURE = %eF10e00% PSI %e/e30Xe "MOST PROBABL 
CE PRESSURE = %oF10.00¢" PSI %e/e30Xe*FORM COEF. (K) =%sF 10.57) 

TIMPLS=0.0 
DO 200 I=1-+eNSTA 

200 TIMPLS=TIMPLS + AIMPLS(T) 
TSLAMT=TIME (NCENTR) «OUR (NCENTR) 
TTT=TIMPLS 
TA © PI®TIMPLS/(2.¢TSLAMT) 
LOCI=NCENTR 
WRITE (MelINI7) TaeTSLAMT 

017 FORMAT (20X-" AVERAGE SLAM FORCE = *eF 8.00% *SINE (PI®T/! oF 6e4e") 

C KIPS*,//) 
WRITE (“eIQIR) LOCI 
018 FORMAT (25xX+*MIDPOINT LOCATION IS HALFeSTATION® «14e///) 
TIMPLS = 0.9 
DO 202 I=I-«NSTA 
202 TIMPLS=TIMPLS*AIMPLS (1) SSHAPE (T) 
TB=P1eTIMPLS/(2.*TSLAMT) 
WRITE (Mel101I9) TBeTSLAMT 
019 FORMAT (20Xe'*MODAL AVERAGE SLAM FORCE = *eF 8,06°*SINE (PI*T/* F644 
Cet) KIPS*#//) 
WRITE (MelN18) LOC] 
WRITE (MeIN69) TTT 
1069 FORMAT (*0%+10Xe*TOTAL SLAM TMPULSE = *.FB.00" KIP-SEC*s///) 
IF (103.£0.0) GO TO If0 
CALL PLOTTING ROUTINE FOR TIME HISTORY OF FORCES 
CALL PLOT (TIMEsFORCE eDELTATsNSTAsIPLToXSCL) 
GO TO 100 
999 WRITE (M-1010) 
1010 FORMAT ('0%s* END OF JOR®) 
STOP 
END 


oak 
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10 


SUBROUTINE UGLYDK (NINeNOUT oO XING YINe XOUT es YOUT sESLeESR) 

THIS SUBROUTINE TAKES THE OFFSFTS AT EACH INPUT WATERLINE AND 

USES SPLINE-CUBIC INTERPOLATION TO DETERMINE THE 10 EQUALLY 

SPACED WATESLINE OFFSETS BETWEEN AL AND 1/10 DESIGN DRAFT 

REAL®R HeNeAeSeAF eBRE oe CE DE HALF oe TWOeSIXeRANSSLPSDOSINeS NCOS 

DIMENSION “IN(S) eYIN(S) eXOUT(1)) eYOUTCI 1) eH OS) 69(5) 6A(25) 6S(5) ¢ 
AE (S) «BE (5S) eCE (5) «NE (5S) 

DATA HALF/N.S0N0/ «TWO/2.0000/+SIX/6.0D00/ eRAD/1 . 7453292590=-02/ 

NMI=NIN-] 

NM2=NM1l-]) 

NEQ=NMe 

00 1 N=leNM) 

HIN) =XIN(N@]) =X IN(N) 

DI(NIECYIN(N* 1) -YIN(N)) ZHON) 

IF «&€SL.GT.90.0) GO TO 2 

NEQ=NE9¢] 

IF (ESR.GE.90-0) GO TO 3 

NEQ=NEOe] 

NSQ=NEQ##?2 

0O & N=] «NSO 

A(N)=0.0 

J=] 

L=) 

IF (ESL.GE.90.0) GO TO & 

A(1)=TWO?H()) 

A(2)=H()) 

SLP=ESL*RAN 

S$(1)=(0(1)-(OSIN(SLP) O0COS(SLP)))@SIX 

J=Je2 

M=LeNEQ 

A(™) =H()) 

L=Me] 

00 S N=leN2 

A(L) =TwO*®(H(N) +H(No)])) 

IF (NeLTeN42) A(L*1) =H(N®1]) 

S( J) =(O0(N*1) -D0(N)) SIX 

J=Je] 

L=Le¢NEQe] 

IF (ESR.GE.90.0) GO TO 7 

A(L-1) =H(NM)]) 

A(L) =-TWOPH(NM)) 

L=L-NEQ 

A(L)=-H(NM)1) 

SLP=ESRERAD 

S(J)=(0(NM])*(OSIN(SLP) ZOCOS(SLP))) @SIX 

CALL SIMO (AeS eNEQeKERROY) 

IF (ESL.LT.90.) GO TO 9 

DO 9 N=l]eN2 

M=NM2-Ne2 

S(M) =S(M=-)]) 

§(1)=0.0 

1F (E€SR8.GT.90.0) SI(NIN) =0-0 

00 10 N=leNM) 

AE (NI =H(SIN@LIKS(INDIS(ISIXPHI(N) ) 

BE (N) =HALF &S(N) 

CE (N) ZO (N) -H(N) @(TWOPS(N) *S(NO1)) ZSIX 

DE (ND SYINON) 

J=)] 


2 





14 
16 


1S 


12 


11 


0O 11 N=lseNOUT 

IF (XOUTI(N) .GE.XIN(2)) GO TO 14 
J=1 

GO TO 15 

IF (XOUT(ON) .LT.XIN(NM})) GO TO 16 
J=NM) 

GO 70 15 

IF (XOUT(N) .GE.XIN(J*})) GO TO 12 
HI) =XOUT(N) -XIN(JU) 

H2=H)] eee 

H3=H)] *He 

YOUT (N) =AE (J) *H3°BE (J) @H] CE (J) 4H) DE (J) 
GO TO }} 

J=Jeo) 

GO TO 16 

CONTINUE 

RE TURN 

END 
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Cc 
C 
Cc 
1002 
C 
C 
C 
19 
20 
C 
. 
C 
3000 


SUBROUTINE SPO0 (XeVKKe89HeSeBBeHHe LAOP SISTA) 
SUBROUTINE SPO00 (XeVKK eBeH eS eRBeHHe JAOPLISTA) 
THIS PROGRAM CALCULATES THE FOR™ COEF. K DEVELOPEN BY OCHT TO SE 
USED IN PREDICTING THF SLAM CHARACTERISTICS OF A SHIP FROM THE 
SHIPS LINES, REFERENCE SNAME TRANSACTIONS 1973 
DIMENSION KCL 1) oe XMOL1) ee YM C12) e-2Z4C11) eROOT(S) o-W(G) o 
XNCL Le ANCL 1) e XL C11) oX2011)9621 0611) 022611) 

oWW (6) eROOTR(S) ¢ROIOTI(S) «COL (6) 
DATA PI/3.1415926/+6FPS/0-000S7 eN/BS MISS 
IF (CTAOP.EQ.1) WRITE (4.1001) Bete BSe4H 
IF (TAOP.FO.1) WRITE (MeLOLI)D (XCT) eFelel)) 
IF (IAOP.FOQ.1) WRITF (421002) ISTA 
FORMAT (28y%4*% HALF<STATION NUMRER *eI3e7) 


CALCULATE AREA ANO MOMENT OF INERTIA BY SIMPSON'#S QULE. 


DO 10 [=2eRe2 

A=ACG,FX(T) +2. exX( lol) 
A=2e*®S/3e*#(ACX(1LIOX(11)¢4e9%X(10)) 

TI=0.0 

AIT=0.0 

DO 20 T=2s8s2 

AI=AlIel. 

TI=ZTI OG. @x(L)*# (ALPS) #92 © 2.8X(141) eC CAT O1.) #5) 8 #2 
AI=AlIel. 

TIS 2c PSS F(T Ge MX (10) 919. 8S) PP 20K (11) %(10.9S) #02) 


DETERMINE A) eA3eAS, AND Ue 


OL=H/B 

IF (ABS(DL-1.) LE. EPS) OL=OL*EPS 
A=A°%2,/P] 

C2=TI*8./P] 

TK1=8942°(1.°0DL*DL**#2) 

TK2=B*()--DL) 

TK 3=-3.°B*(1.°0OL) 

TKG=(ATKI]I STK? 

TKS=TK3/7TKP 

TK6H=2./TK2 

D1i=N0L*%*%2*4.70L-7. 

O2=-OL%*#3*NLe%2? -3.*%DL+3- 

03=-DL+S.- 

04=D0L**3*0L*®*%2+3.°0L*3. 
05=-0L*%*2-9.*DL-le. 

06=3-.*NL*7. 

07=-1.75°DL %846-0.590L 99 3-0L992-1.590LK 1.25 
W(1)=-4,7%TK6 

W(2)=2.°ReTK69D 3-4. 9TKS—-66 
W(3)=B%(2.%(D66T*S29D3) eRPTKO9D] “4.9 TKS 
W(4)=B8(2.°TK49D 3°R8 (N56 TKS8D1 ©0.59R*TKE%D2) ) 
w (5) =Bee2%(TKG2D)] 9.5920 (3.°04+TKS*02)) 

W (6) =0.59R283% (T4820 2+R907) -Ce 

IF (TAOP.EN.1) WIITFE (4.3000) (W(T) 6 T=106) 
FORMAT (° W= 46(F£15.543X%) 9/7) 

MM=S 

00 30 1=1+4 

J=7-I 

wwt(T)=wW(J) 
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30 


40 


37 
38 


39 


3001 


3002 


CONTINUE 

CALL BAIRS (WwWeCOL eMMeROOTRsROOTI «KK) 
ERROR=0.00000) 

NR=0 

00 40 1=]-5 

ROOT(1)=0.0 

IF (ABS(ROOTI(1)).GT.FRROR) GO TO 40 
IF (ABS(ROOTRI(1I)).LT-FRROR) GO TO 40 
NR=NRe } 

ROOT (NR) =ROOTR(I) 

CONT INUE 

IF (NR.FQ.1) GO TO 39 

00 38 J=1e3 

KB=J¢] 

DO 37 I=KB.NR 

IF (ROOT(J).GT.ROOT(I)) GO TO 37 
R=ROOT (JY) 

ROOT (J) =ROOT(I) 

ROOT(1I)=R 

CONTINUE 

CONT INUE 

IF (NR.FQ.4%) GO TO 39 

IF (ROOT(4).GT.ROOT(S)) GO TO 39 


R=ROOT (4) 

ROOT (4) =ROOT(S) 

ROOT (S)=R 

CONTINUE 

IF (TAOP.EQ.1) WRITE(M,.3001) (ROOT(I) e1=1 eS) eNRekr 

FORMAT (® ROOT = %eSlF190eSe3X) o/e* NR = Ce150% KK = %e15e/) 


IF (KK.GE.}) GO TO 99R 

DO 25 [=]+.Ne 

U=ROOT (I) 

ALF A=2.°eU/9 

AS=TKG/SU4CTKS*TK68U 

Al=(1.-OL)/ALFA-AS 

A3=().°DL)ISALFA-}. 

IF (ITAOP.FQ.1) WRITE (.3002) AleA3,A5 

FORMAT (° Al = %eF1l9.Se3XetA3 = *eF TOS e 3Xe*AS = %0F10.547) 

IF (A3 «LT. 92) GO TOD AS 

IF (ABSCAI-FPS) LE. 1. «AND. ABS(A3-EPS) .LE. Ie AND. 
ABS (AS-EPS) etEe Je) 69 TO 26 


25 CONTINUE 


AAO 


26 


3005 


GO TO 998 
COMPUTE K-VALUF 


TK=-3.59893 © 2.419RB8%AI - 0.8728558%43 * 9,62395%AS 
VKSEXP(TK) 

VKK=VK 

IF (ITAOP.GT.0) WRITE (Me2000) AleA3BeADeUl VK 

IF (IAOP.GT.0) WRITE (Me2004%) 

A=A®PI/4,. 

IF (IAOP.E9.1) WRITE (423005) A 

FORMAT (1X6 HALF<SECTION AREA OF 1/10 DRAFT = *yFI0.Se/) 
IF (VK.LFE.L0.19) GO TO 109 

VKK= 0.027°5*B/s4 

IF (IAOP.E9.1) weITe (Me 3004) VKK 
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3004 FORMAT (1Xe* THE REGRESSION ANALYSIS * NOT USED ---=te/olXeo* K = # 
CoF10.603X%0e" WHERE K = 0,027 *B°B/A?,/) 
109 CONTINUE 
GO TO 100 
998A wRITE (Me2NN2) 
100 CONTINUE 
RETURN 
1001 FORMAT (1H1.S8X.12eHe#* SPOO 848 772BX.72HTHIS POOGRAM COMPUTES T 
2HE CONFORMAL MAPPING COFFFICIENTS AND THE K-VALUE/2B8X. 
2 63HNEEDED TO COMPUTE SILAMMING PRESSJRE BY THE REGRESSION EQUATION 
2e SLS2RX ce LGHHALF SREANTH = eF 10.5eSH FEET e 3X0 
2 2SHHEIGHT ABNVE THE BOTTOM = «F10.5e5H FEETe/28X614HHALF BEAM 
2= oF 10.6595SH FEF Te 21 Xe THNRAFT = oF10.595H FOFT) 
1011 FORMAT (7/7 2RXeIIHTHE VALUES KX eX ecccoce X ARF .FI1O0.S/40Xel1H0> 
2 BKelH le Xe 2HLIeSKe GIF IL0O6S/59X) 0F192595SH FEET) 
2000 FORMAT (/S//27AX e7TGHTHE COEFFICIFNTS OF THE TRANSFORMATION OF THE C 


2IRCLE IN THE ZFTA=PLANE ARESS2ZRXeSHA = 0E15-8e7XeSHA = 615.8% 
2 7X%.SHA = eF YG eAS2ASIK ce PH] oe 25K ea LHI e2P6Xe LHSSSSCAKe LQHTHE SCALF. FACTO 
OR 18 e41 Xs] GHTHE KeVALUE IS//28XeGHU = 0E15,8s36X%,4HK = .F15.A) 


2002 FORMAT (///78%e6K6HA SOLUTION FOR U AAS NOT BEEN FOUND, OLFASE CH 
2ECK THE INPUT DATA.//) 


2004 FORMAT (/36%«¢S7HWHERE K = EX® (= 3.599 * 2.420 A - 0.873 A ¢+ 9 
2-624 A ) ZO9XelLH1Le LOXe LL 43910X%e1HS) 
END 
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SUBROUTINE SLAMPR (RM oRVeHF eHAcALWLE e ISTACEPRES e VK eSLAMS ePRORPR, 
DVEL eVTHRSe TSTORMeHSTGePRSLAMePVEL eSTASP) 

THIS SUBROUTINE CALCULATES THE EXTREYE PRESSURE AT EACH STATION 

WHERE THE PROPABILITY OF SLAM IMPACT IS GREATER TUAN } 

DIMENSION R4(17) -RV(17) eEPRES (17) eVK (17) ePROBPR(17) eDVEL (17) ¢ 


CPRSLAM(17) «SLAMS (17) ePVEL (17) 0X (17) 


DATA ALPHA 0&2 eP1/0.01632.17403.1415926/ 
DATA MeN/5eR/ 

RSTARS=VTHRS#e2 

AL=0.983°ALWL 

T=TSTORM 

NSTA=ISTA 

NO 10 T=1]eNSTA 

H= ((HA=HF) 740.) °FLOAT(]) oHF 
DRM=2.%( (AL PRAM(T) 90.83) e%2) 
DRV=2.°GPAL% (RV(IT) 20.83) 992 
PRSLAM(T) =© x08 (-HOH/DRM-RSTARS/SORV) 

SLAMS (1) =(3400.%T/(2.%PJ)) ®SORT(DRV/IRM) PPRSLAM(T) 
IF (SLAMS(T).GE.1.) GO TO 2) 

ISTA=I-) 

NSTA=ISTA «+ ] 

WRITE (Mel19N2) NSTAcTSTORM. SLAMS (1) eASIG 


1002 FORMAT (*0%.* THF NUMRFR OF SLAMS AT ALF STATION *el3e* NURING A 


i1 


10 


CSTORM OF DURATION oF 3.00" HOURS = teoF7e4e/e* FOS A SIGNIFICANT W 
CAVE HEIGHT OF %eFS5.00* FEET* eo //) 


GO TO 100 

EX=]-/SLAMS (1) 

EPRES(II=AVK (I) @(RSTARSHNAVEALOG (1 e—- (1) e-ALPHA) @OFX) ) 
PROROR (TI=vKE(TI*®CRSTARS * DRV®ALOG(SLAMS(I1))) 

OVEL (1) =SORT(EPRES( I) /VK(T)) 

PVEL (1) =SOR8T(PRORPR(I) sVK(1)) 

CONTINUE 

ISTA=NSTA 

WRITE (Me]NN]) ISTAeSLAMS(ISTA) eTSTORMSHSIG 


1001 FORMAT (°1%.* ADDITIONAL STATIONS SHOULD BE CONSINEPRED SINCE AT HA 


100 


CLF STATION *%el3e%e* THF NUMBER OF SLAMS = #,F5e.)6* FOP A STORM OF 
COURATION *.F5.00% HOURSt#e/e* THE SIGNIFICANT WAVE HEIGHT = %eF 5.0% 
C* FEET*.//) 


RETURN 

CONTINUE 

NCENTR=ISTA/S2¢]) 

HSTASP=STASP/2. 
XCENTR=FLOAT I (NCENTR) @HSTASP-HSTASP/2 6 
DO 1 J=leNCENTR 

X(J) =XCENTR=- (FLOAT (J) PHS TASPHHSTASP/2.~) 
CONTINUE 

J=NCENTR*] 

00 2 I=JeNSTA 

X(T) =X(T-1) *HSTASP 

CONTINUE 

DO S IT=leISTA 
A=-112.eHSIG*x (1) / (AL #2) 

FPRES(IT) =EPRPES(I)2®(10.204) 
PRORPR( IT) =PROBPR(I) (10.794) 

CONTINUE 

RE TURN 

ENO 


D7 








SUBROUTINE RELMV (RM RV) 
THIS SUBROUTINE TAKES THE INPUTED RELATIVE MOTIONS AND VELOCITIES 
AT STATIONS 1 & &© AND ASSUMES LINEAR VARIATION TO YETERMINE VALUES 


aAaAO 


10 


20 


AT THE 17 HALF=-STATIONS 
DIMENSION PM(17).RV(17) 
SLOPE=(RM(8) =-RM(?))/6, 
B=RM(2)-SLOPE?2, 

DO 10 T=1l-17 

RM( IT) =SLOPECFLOAT(I) +B 
CONTINUE 

SLOPE = (RV(8)-RV(2))764. 
B=RV(2)-SLOPE?, 

00 20 J[=1-17 
RV(I)SSLOPFOFLOAT(I) «© AR 
CONTINUE 

RETURN 

END 
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SUBROUTINE FORCES (FORCE oNSTAso TIME eDELTAT ePRESoXX,ALBPeOVEL o VAFTS» 
C STASP eVYePVEL + I0O4) 
THIS SUBROUTINE DOES THE SPATTAL AND TEMPORAL INTFGRATIONS AT EACH 
HALFeSTATION OF THE PRESSURES TO FIND THE SLAM FORCE TIME HISTORY 
OIMENSION FORCE (17020) ePRESCI7ol lL) oe XXC17ell) oe TIME (17) eOVEL(17) © 
c YY(17) ePVEL (17) »© (10) oF OR(17.20) 
TSLAM=0.18SQRT (ALBP/570.) 
OELTAT=TSLAM/19. 
NOELT=INT (STASP/(2.*VAFTS*®DOELTAT) ) 
IF (NOELT.LT.1) NOELT=) 
T2STASP/(VAFTS#2,) 
NCENTR=NSTA/2¢) ions 
TIME (1)=0.0 . aa 
DO 4 I=2eNCENTR 
TIME CI) =TIME (Il-1l) ¢T 
CONTINUE 
IZNCENTR*] 
DO S J=IeNSTA 
TIME (J) =TIME (Jel) oT 
CONT INUE 
00 20 I=leNSTA 
IF (106.£Q.1) GO TO 50 
OELY=DELTAT®0OVEL (1) 
G60 TO 5) 
SQ OELY=OELTAT®PVEL (I) 
Si FORCE (1e+1)=0.0 
00 21 J=2e5 
21 FORCE (LeU) = (FLOAT (Jol) s5e) (PRES (Tel) OXXC1o1) © PRES(T 52) @(XX(1e2) 
CoXX(101))) 
DO 22 J=6910 . 
22 FORCE (Le J) =(FLOAT(L1<“U) 750) @(CPRES (Tol) OXX(T41) © PRESCTe2)* (Xx (Te2 
Cy-XxXC1e1))) 
00 28 J=11-+20 
28 FORCE (1-J)=0.0 
00 23 J=3e)11 : , . 
F(1)=0.0 "Ys 
OELX= (XX (Te J) exx(TeJol)) @PRES (I J) ” 
DO 24 JJ=2e5 : -- 
26 F (JU) = (FLOAT (JU=1) 756.) @NELX a - 


00 25 JJ=6010 | | fal ee: 


2S F (QU) = (FLOAT (C1 1-JU) 75.) @OELX EUG Ae ee ae 
00 30 L=1-20 , 0G 
MINT (YY (J-1)/ (OVEL (1) @FLOAT (L) *OELTAT)) 
INCeL 
IF (MeLlE-1) GO TO 3) 
30 CONTINUE 
GO TO 20 
31 M=INC 
MM=M¢9 
IF (MmM.GT.20) GO TO 100 
DO 26 JJ=1-10 
K=MeJJ=-])] 
26 FORCE (I eK) EFORCE (I oK) ¢ F (JU) 
G60 TO 23 
100 IF (M.GE.20) GO TO 20 
mM=20 
K=0 
00 27 JJ=M.MM 
K=Ke} 


a9 





27 FORCE (1 +JJ) =FORCE (16 JU) oF (K) 
23 CONTINUE 


29 
20 


80 
70 
60 


00 29 J=1.20 

FORCE (I -JI=(STASP®FORCE (1¢J)/1000.) ©144./FLOAT (NDELT) 
FOR (Ie J) =FORCE (I eJ) 

CONTINUE 

IF (NDELT.LE.1) RETURN 

00 60 T=l!sNSTA 

NEXT=)} 

LAST=0 

00 70 M=2.NDELT 

NEXT=NEXT+} 

LAST=LAST+} 

DO 80 J=NEXTe20 

JJZI“LAST 

FORCE (I oJ) =FORCE (Ie J) FOR (Te JU) 
CONTINUE 

CONT INUE 

CONTINUE 

RETURN 

ENO 


100 





aAaAaAa 


¥2 
13 


14 
10 


SUBROUTINE PRESS (EPRESsNSTAsPRES o XX 0 VV eS eXG) 


THIS SUBROUTINE DETERMINES THE SPATIAL DISTRIBUTION OF THE VERTICAL 
COMPONENT OF THE PRESSURE AS A FUNCTION OF DISTANCE FROM THE CENTESLINE 


DIMENSION EPRES(17) ePRESCI Tol 1) oXX(L7oLL) eVVC11) eCTHETACIO) eXC11) 


CoXG(17e1)) 


DO 10 T=leNSTA 

PRES (1Te1) =EPRES(T) 
XO1L)=XX(1e)) 

DO l2 v=) <10 

DEL TAX=XX(TeJel)—-XX( Te) 
HYP=SORT(S®e2eNELTAx##?) 
XC Jel) =X J) +HYP 

CTHE TAC J) =9ELTAX/SHYP 

00 13 J=le1)} 

XGCTeJP =X) 
DELTAP=EPRES(1I)/10. 

DO 14 J=2e)) 
A=FLOAT(J-1) 

PRES (IT eJI =H=(EPRES(I) -A®NELTA®) SCTHETA(J=)) 
CONTINUE 

RETURN 

END 


coy 








20 


aAaAQD 


)) 
le 
13 


10 


SUBROUTINE KEFFEC (VKeAVGK) 

THIS SUBROUTINE CALCULATES THE AVERAGE OF THE FIRST 11 HALF- 
STATIONS CALLED THE AVEPAGE FORM COEF. 

DIMENSION vK(17) 

AVGK=0,0 

00 20 J=1lel) 

AVGK=AVGK «VK (J) 

AVGK=AVGK/}]. 

RETURN 

ENO 


SUBROUTINE APPROX (NSTAeFORCE oF MAX eAIMPLSe Ae DUR sDFLTAT) 
THIS SUBROUTINE CALCULATES FOR EACH 4ALF-STATION THE IMPULSE RY 
INTEGRATING THF FORCE TIME AISTORY AND COMPUTES A HALF-SINE 
APPROXIMATION WITH THE S4ME IMPULSE AND DUSATION 

DIMENSION FORCE (17¢20) e0FMAX(17) eATMOLS(17) 64017) eNUR(17) eV (20) 
DATA PI73.141592654/ 

oOo 10 Iz]eNSTA 

00 11 J=10.70 

IF (FORCE (TeJ).GE.0.05) GO TO 11 

DUR (IT) =FLOAT (J-1) *DELTAT 

NPTS=J 

GQ TO le 

CONTINUE 

OUR CL) S20. °DEL TAT 

00 13 J=1.20 

Y (J) =FORCE (1 eJ) 

CONTINUE 

CALL MAXFOQ (Y.AMAX) 

FMAX (1) =AMAX 

CALL OSF (NELTAT+YsYe2N) 

AIMPLS(T)=Y(20) 

A(I) = PIeAIMPLS(I) 702. 290UR(T)) 

CONTINUE 

RETURN 

END 
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Aa 


AAA 


10 


11 


20 
10 


SUBROUTINE MAXFOR (YeAMAX) 
THIS SUBROUTINE DETERMINES THE MAXIMJM FORCE FROM THE FORCE TIME 
HISTORY AT EACH HALF=<STATION USING A PARABOLIC FIT THROUGH THE 
HIGHEST 3 POINTS, 
DIMENSION Y(20) 
FLAST = Y(2) 
00 10 T=348 
N=I 
IF (YCT) eLT-FLAST) GO TO 1) 
FLAST = YCT) 
CONTINUE 
AMAX = 1.0£F10 
AA=Y (Ne-2) 
BB=Y (N-]J) 
CC=Y(N) 
B=0.5% (4. *AR-3, SAA-CC) 
C=AA 
=BR-AA-B 
K=—-B/ (PFA) 
AMAX=A®X8OP «¢ ROK « C 
RETURN 
END 


SUBROUTINE PLOT (TIME FORCES eDEL TAT eNSTAcIOL T+ XSCl ) 

THIS SUBROUTINE USES THE MeievTe SeCleFe PYICTR PLOTTING ROUTINE 
- IF THE PSOGRAM IS TO RE USED ON 4NOTHER MACHINE. THIS IS THE 
ONLY SUBROUTINE THAT REQUIRES CHANGING ALONG WITH THE RELATED 
CARNS IN MAIN 

INTEGFRe2 TXPLOT(17) «xt a3(40) 

DIMENSION PL TR (34020) e TIME (17) eFORCES(17e20) oe XSCL (G) 

DATA IXPLOTZ1I Re 1D e202] 022 0230240250280 279 28029030031 0320330347 
00 10 J=1eNSTA 

T=TIME (T) 

00 20 J=1-20 

PLTR(IeJ) H=FORCES (I « 3) 

Tl=1I¢17 

PLTR(IIT«J) =T 

T=T*DELTAT 

CONTINUE 

CALL PICTR (PLTRe34eXLABOXSCL eNSTA 20 emlemlLeSdeIPLTeOele I XPLOT) 
PAUSE 

RE TURN 

END 
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SUBROUTINE BAIRS (ACOF «As NORDsROOTReROOTI +s IER) 
REAL ACOF (1) +A(1) eROOTR(1) -RONTI(}) 
DATA TOL2/) .0F -3/sTOLX/S- OF -6/ 
DATA TOL+ITMIN,ITMAX/1.05-9910+100/ 
DATA TOL3/1.NF-10/ 
N=NORO 
IF (NeLT.LO)N=-N 
O=ACOF (Ne) ) 
00M=-100.0 
IR=0 
Ce-----BASIC ERROR CHECKING 
IF (ND Lele? 
1 JR=] 
GO TO 200 
IF (N-99)4%.4,3 
TR=2 
GO TO 200 
IF (0)60536 
TR=4 
GO TO 200 
00 10 T=leN 
0 A(T) =ACOF(1)70 
A(Nol]l)=1.0 
REMOVE TRIVIAL ZEROES FIRST 
J=2 
00 14 J=leN 
IF CACT)) 12061112 
}) ROOTR(J-1)=9.0 
14 J=Je} 
Ceo---- FIND A REAL ROOT IF REDYCED POLYNOM IAL IS O00 
l2 LEF T=N-Je) 
IF (LEFT) 100-90.13 
13 NN=Ne}] 
I J=J-1 
Ce----- ESTIMATE THE MODULUS OF THE LARGEST ROOT 
00 151 IT=IUeNN 
IF CACTI) .NE.O0.0) 
IDOOM=AMAX) (AL OG (ABS (4(T) ACT I) )) FLOAT (IT -14*]) eD0M) 
15] CONTINUE 
DOM=ExXP(D0™) 
QJ=)].0/D0M/ND0M 
PJ=-2.0/00™ 
G SCALE THE POLYNOMTAL SO LARGEST ROOT IS NEAR 1.0 
00 9 J=TIJeNn 
00 9 K=I Je! 
9 A(K)=A(K) /0N0M 
IF (LEFT.NE.LEFT/S2%2)60 T0 15 
300 x0=0.0 


WN 


om Wa & 


PO=A(Jq-}) 
Cee----DETERMINE IF REAL ROOT IS POSITIVE OR NEGATIVE 
X1=1.0 
IF (P0.GE-0.0)X1==-1.0 
IP=] 
Coren EXPAND RANGE UNTIL ZERO HAS BEEN SRACKETED 
c 
Coreoe-- THIS IS AN INTERNAL SUBROUTINE TO EVALUATE A REDUCED POLYNOMT AL 
fc 


400 ICNT=Ne2 


104 








4} 


43 


14] 


30 


P1=0.0 
PI=P1l®X1*A(ICNT=-1) 
ICNT=ICNT-} 

IF CICNT=-J) 4022401 640) 
GO TO (320.330).IP 


IF (PO®P1) 3254360432) 
x0=X)} 

X2=10.00X] 

GO TO 400 


~-ANSWER LIES BETWEEN X09 AND X) 


IP=2 

X2=X] 

X1=(X0¢X2) 290.5 

IF CABS (X2-X0)-ABS(X]) *TOLX) ) 360.360.2400 
IF (PO*®P1) 331 «360332 

X2=X]1 

GO TO 376 

x0=X] 

GO TO 326 


--NOW DIVIDE OUT THE ZERO 


ROOTR(J-1)=x14N0M 

]=N 
ACTI=A(T) ox] °A(T¢)) 
l=!I-] 

IF (1-J)380-.370.370 
J=J¢] 

Conn. Uc 

NO BO J=IJeNeo2 

ITRY=0 

IF (JeLTLN)GO TO 4] 

P=A(N) 

Q=A(N-1]) 

GO TO 50 

IF (NORD.GT.9)IGO TO 43 
TRY PREVIOUS ROOT FOR STARTING POINT 

P=PYF®ROOTRIJ) 

Q=O4* (ROOT2 (J) ®ROOTR (YJ) + 2OOTI (YU) FROOTI (UY) ) 

GO TO 141 

P=0.0 

O=1TRY 

IF CITRY.GT.1)GO TO 47 
ITRY=ITRY >} 

DO 45) ITES=) +I TMAX 

BO=0.0 

B1=0.0 

C0=0.0 

C1=0.0 

C2=0.0 


--COMPUTE COLUMN OF BtS AND C#S 


I=Ne-] 
OD=A( 1) -P #8] -Q#R0 
BO=B) 

B1=00 

l=lI-] 

IF (1-JU#2) 40240430 
ODO=Bl-P#C2-0*C) 
CO=C) 
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Cl=Ce 
Cc2=00 
GO TO 20 
Conwnwnm--SPECIAL FORM FOR LAST COEFFICIENT 
40 C2=-P"#C)]-0°C0 
OD=C!}"C}-C2eco0 
DP=B0"C}-81*C0 
DO=R1"*C)])-RNeC2 
IF (A8S (DD) .LE.-ABS(TOL2*DP) .OR. ABS (DD) -LE~ABS(TOL2*DQ))G0 TO 43 
DP=DP/DD 
DO=DQ/DD 
P=P+DP 
Q=0-D0 
ERR=DP#DP+N0eDQ 
CRIT=P9P+OQ80 
IF (CRITLLE.TOL2®ERR)IGO TO 43 
CRIT=ERR/CRIT 
IF CITER.GT.3.4ND.CRIT.LYe¢TOL3)GO TO 50 
Conmn----9N0 AT LEAST ITMIN ITERATIONS 
IF CITERP<$ITTYIN)D 454645244 
Coeenwnm--ITF ERROR 15 DECREASINGe KEEP GOING 
44 IF (CCRIT.LGE.CRITO)GO TO 46 
IF (CRIT.LLT.LTOL3IGO TO SN 
45 ERRO=ERR 
CRITO=CRIT 
465) CONTINUE 
Cemme-=-NROP OUT wHEN ERR STARTS TO INCREASE OR ITER IS GREATER THAN ITMAY 
46 IF (CCRITO-TOLI59050e43 
Cowwnnn]F ACCURACY IS INSUFFICIENT. GIVE FRROR OUTPUT A*ID PPOCEFD 
47 IR=3 
WRITE (SeSONOOIEPROCERReCRITOSCRIT 
S000 FORMAT(1P.19G)]2.4) 
Ss0 pP=—_P/2?.,0 
DISC=PPePP-Q 
IF (DISC) 5245454 
Se DISC=SORT(-DISC) 
ROOTR(J-I) =PPeNOM 
ROOTR( J) =P9°eD0™ 
ROOTI (J-1) =-DISC?D9M 
ROOTI (J) =DISC#NOM 
GO TO 56 
54 DISC=SQPT (DISC) 
ROOTR(J-I)=(PP=-DISC) *DOM 
ROOTR( J) = (PP *+N1TSC) #DOM 
ROOTI(J-1) =0.0 
ROOTI(J)=0.0 


C-----= SYNTHETIC DIVISION WITH QUADRATIC FACTORS 
56 I=N 
JL=J 
A(N) =A(N) =P 
60 T=I-] 


IF (IJ) 80+70+70 
70 A(ITI=A(I) -PPA( TO} )-“OQ9A( TP?) 
GO TO 69 
80 CONT INUE 
IF (UL-N) 902-100-100 
90 ROOTR(N) =-<A(N) 
ROOTI(N) =0.0 
100 CONTINUE 
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200 IER=IR 
RETURN 
END 
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APPENDIX C ~- PROGRAM USER'S MANUAL 


This appendix will discuss the calculations performed in 


the various subroutines, the output format options, and the 


input data deck. 


DESCRIPTION OF PROGRAM 

MAIN ~- The Main program does some simple calculations, 
contains formats for I/O, determines from the I/O options 
specified what calculations are to be performed, and calls 


the appropriate subroutines to perform the calculations. 


SUBROUTINE UGLYDK - This subroutine takes the five 
inputed offsets and uses a Spline-cubic curve fitting routine 
to interpolate ten offsets in equal increments of T'/10 from 
the base line to T'. The subroutine is called once for each 
set of inputed offsets. The subroutine requires solution to 
a set of simultaneous linear equations which is done by the 
IBM SSP/SIMQ subroutine which should be available to most 


users. 


SUBROUTINE SPOO - This subroutine performs the calculations 
necessary to find the form coefficient k. For values above 
k = 0.19 equation (6a) is used to determine k. This subroutine 
performs the mapping necessary to determine the regression 
analysis form coefficients. This subroutine is called once 
for each set of inputed offsets. The subroutine requires the 
roots of a fifth degree polynomial which is accomplished in 


the BAIRS subroutine. 
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SUBROUTINE SLAMPR - This subroutine takes the relative 
motions, relative velocities, form coefficients, Significant 
wave height, and drafts and calculates the following: 

1. number of slams at each half-station, 

2. longitudinal extent of slan, 

3. probability of slam at each half-station, 

4. most probable extreme pressure at each half-station, 

9. extreme pressure for design at each half-station, 

6. and the most probable and extreme velocities of 

impact at each half-station. 
This subroutine uses the various equations developed in Section 


IV to determine the above information. 


SUBROUTINE RELMV - This subroutine takes the non-dimen- 
sional relative motions and relative velocities inputted at 
station 1 and 4 from the tables and uses linear interpolation 
to find the values for the 17 half-stations. Seventeen half- 
stations were considered since slam damage has been reported 
as far as 40% aft of the forward perpendicular. Linear inter- 
polation is valid within the linear theory used in seakeeping 


programs. 


SUBROUTINE FORCES - This subroutine performs the temporal 
integration of pressure determined by SLAMPR and PRESS sub- 
routines. The output is the pressure time history for each 


half-station in the slam region. 
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SUBROUTINE PRESS - This subroutine takes the pressure at 
the keel from SLAMPR and calculates the girth distribution of 
vertical pressure assuming linear variation from the flat to 
T' and resolves the vertical component at each elevation using 
the hull offsets to calculate the angle 8 that the normal 


makes with the vertical (Figures 8a & 8b). 


SUBROUTINE KEFFEC - This subroutine is not used in any 
of the calculations but determines an average hull form coeff- 
icient. This value is a measure of a ship's overall suscep- 
tibility to slam forces. For two similar ships the one with 
the larger average k will have larger slam forces and impulses. 
The average is done over the forward eleven half-stations 


since slamming rarely occurs aft of this point. 


SUBROUTINE APPROX - This Subroutine uses Simpson's 
integration to find the slam impulse by integrating the 
pressure time history from FORCES at each half-station. This 
Subroutine calls the IBM SSP/QSF subroutine to perform the 
integrations. Once the impulse is known, the amplitude of 
an equivalent half-sine pulse with equal duration can be 
found. This sine approximation is determined to input into 


the S.H.V.R.S. dynamic analysis program. 


SUBROUTINE MAXFOR - This subroutine fits a parabola 
through the three highest force points at each half-station and 


the first derivative is set equal to zero to find the maximum 


force at each half-station. 
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SUBROUTINE PLOT - This routine uses the M.I.T. Joint 
Computer Facility PICTR plotting routine to plot the force 
time history for each half-station on one graph. Figures 9 
Bente are computer plots using this routine. This is the only 


Subroutine that would require changing if running the program 


on another computer. 


SeTPUL OPTIONS 
There are several output options available from the 
program depending on the information desired. Options are 
specified in the data deck by assigning integer values to 
integer variables. The following options are available: 
fiero = | This Option prints information from the SPOO 
subroutine to aid in determining problems in this 
routine. See Figure 22 for the information printed. 
T0l1 = @ Delete the above output. This is the 
normal option. 
yeeros = 1 This option allows a plot of the pressure 
time history. See Figures 13 to 16 for typical plots. 
103 = @ Delete the plots. 
3. 104 = 1 This outputs a table of inputed offsets and 
the ten calculated offsets for each half-station. It 
is desirable to specify this option the first time a 
set of lines are read so as to check the data. See 
Figure 23 for an example output. 


T04 = @ Delete table of offsets. 
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I06 = 1 This option specifies that the calculations 
and output values will use most probable extreme 
pressure values. 

I06 = 9 This option specifies that the program uses 
the extreme pressure for design. 

[07 = 1 The force time history is outputed in 
tabular format. See Figure 24 for an example of 

the output information. 

I07 = @ Delete table of forces. 

I01f = 9 Write table of information for each half- 
Station in the slam region. See Figure 25 for the 
information available in this option. 

TO1f = 1 Write a single table of information for 
all the stations in the slam region. Tables 7 to 10 


were prepared using this option. 


I019 = 2 Delete both of the above tables. 


Any combination of the above options may be Specified to get 


the information desired. Since the main cost of running the 


program is the output cost, it was felt that being able to 


eliminate unnecessary output information was necessary. 


Several input options are available: 


Ile 


105 = 1 This allows the inputting of all the relative 


motions and velocities from seakeeping results if 


available. 
I05 = @ Input relative motions and velocities for 
stations 1 & 4 only. (Non-dimensional values from 


Tables [15]). 
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Time of storm may be specified if desired. If 

TSTORM is specified, it will be used in the calcul- 
ations. If set equal to zero, the TSTORM values will 
be obtained from Figure 4. 

The traveling velocity may be specified if desired. 
The value specified should be for a 520 foot ship 
and the program will Froude scale this figure. If 
the traveling velocity (VAFT) is set equal to zero, 
the traveling velocity will be calculated by the pro- 
gram as discussed in Section IV. 

The Froude number or the speed in knots can be 


specified. Both need not be inputed. 


DATA DECK 


In t 


be discus 


CARD SET 


CARD 


CARD 


his section the input format for the data deck will 
sed. 
# 1 
1 = Format (3 I 5) - I01, 102, 104 
The values for I0l1 and I04 are discussed above. 
102 is always set equal to 1 and indicates that 
ships lines follow. 
2? - Format (20 A 4) - Title and information on ship. 


This card contains any information desired by the 
user to designate the ship being run. This data is 
printed in the output as a heading only and is not 


used by the program. 


ie 





CARD 


oe FOrmat (SeFs10.5)m@— SbaehH ALBP, ALWL, STASP 
BB = Half Beam 
HH = Design Draft 


ALBP 


Length Between Perpendiculars 


ALWL 


Length of Load Waterline 


STASP = Station Spacing (not half-stations) 


CARD SET #2 


This 


card set contains the ships lines necessary for 


slam calculations. 


CARD 


CARD 


CARD 


lien Hormac. (5 fF 10,5) 

This card contains the waterlines at which the 
offsets were taken. The first value must be 
the 0.0 foot waterline. 

Jae ORMaces( 5. F 410, 5) 

The offsets for the five waterlines in the same 
order as CARD 1 for the first half-station. 
CARD 2 is repeated for as many half-Stations as 
desired up to a total of 17 half-stations. 

3 - BLANK CARD 

This card indicates that the last set of offsets 


has been inputed. 


CARD SET #3 


This set contains only one card. If no plots are desired 


(103 = Q@), 


this card is left blank. 
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aomiacw (ia, 4 F 10.5) = TPLPA XSCL 


iPGie> +l — This specifies that the plot routine 


scales the plots. (AUTOSCALE) 
bias ~2 = thes specities tovthe plot routine that 


the values in XSCL be used to scale the plots. 


XSCL(1) = 0.0 = Time slam starts 

XSCL(2) = End of time axis 

XSCL(3) = 0.0 

XSCL(4) = Maximum force on vertical axis. 


If IPLT = +1, XSCL need not be specified. This same 


data will 


are read. 


be used to scale all plots until a new set of lines 


CARD SET #4 


Card 


set #4 contains the data necessary for a particular 


set of operating conditions. This set is repeated for as many 


operating 


conditions as desired. 


opp le- Format (6 I 5) saei09, 103, 105, 106, 107, 1019 


CARD 


T09 = 1 always = 1 - indicates a new set of operating 
conditions follow. The other integer variables 

were discussed above. 

2 - Format (20 A 4) 

This card contains information which is printed on 


the output to distinguish this operating condition 


from the others. 
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een 3 > Format (8 F 10.5) - FROUDE, SPEED, SEAS HF, Hae 


CARD 


CARD 


CB, VAFT, VTHR 
FROUDE = Froude number 
SPEED = Ship speed in knots 
SEA = Significant wave height (feet) 


is 


Draft Forward (feet) 

HA = Draft Aft (feet) 

CB = Block coefficient of underwater portion of 
hudivior this draft. 

VAET = Traveling velocity for 520 foot ship (ft/sec) 

yun —~siheeshold velocity for S20mtoot ship (ft/sec) 

4 - (F 10.5) - TSTORM 

TSTORM = Duration of storm in hours 

> - IF I05 = @, Format (4 F 10.5)- RM(2), RV(2), 


RM(8), RV(8) 


RM(2) = relative motion from tables at station l 
RV(2) = relative velocity from tables at station l 
RM(8) = relative motion from tables at station 4 
RV(8) = relative velocity from tables at station 4 


feels —= 1 Format (8 F 10.5) for each card = use 

aS many cards as necessary to specify the non-dimen- 
sional relative motions and velocities for as many 
half-stations as offsets read in Card Set #1. The 
values are specified in pairs beginning at half- 


Seaton I (1.e.,;, RM(1),  RV(E), RM(2);, RV(2)>.....- ) 


ake 





CARD 6 - IF I03 = 9 -- Delete Card 
IF 103 = 1 -- FIRST 40 spaces contain X-axis 
label 


NEXT 40 spaces contain Y-axis 


label 


This card is used by plotting subroutine. 


At the end of the operating conditions for this ship a 


blank card is inserted to indicate the end of 
for this set of lines. The above card sets l 
repeated for as many ships as desired. After 
has been read, another blank card is inserted 
end of the data deck. Therefore the last two 
will be blank cards indicating the end of the 


tions and that the last ship has been read. 


ey 


Ene wcalculations 
to 4 can be 

the last ship 

to indicate the 
cards in the deck 


Operating condi- 
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APPENDIX D - EXAMPLE PROGRAM 
An example program data deck and sample output are 
included to aid in using this program. MARINER case 2 is 


the example run which is the fully loaded case at 15.5 knots. 
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4/7 XEQ SLAMCW 


0 1 0 
MARINFR = = = = = s+ sf ee ef |= = - - AS BUILT 
38.0 29.75 528.0 537613 266 
0.00 1.00 2200 3.00 4.00 
0.5 1.5 1.75 2.05 2.3 
1,00 3.00 Se > &.) 4.5 
1,0 3245 4.3 5215 5.6 
be 0 3.9 Sel 6.2 6.6 
15 4.95 6.3 7.5 8.05 
2.6 6.0 70% 8.8 9.5 
4.0 7.5 9.2 10.6 11-6 
5.4 9.0 10.9 12.4% 13.5 
6.95 11.0 13.1 14.8 14.0 
8.00 12.00 15.00 ie 18.2 
8.58 1S. 17.6 19.7 21.9 
9.17 18.00 20.2 22.3 24.08 
13-6 21. 23.10 24.9 26-5 
RL AN 
-2 0.90 0.5 0.90 S000. 
1 1 0 0 0 1 
2 MARINER = F=H0.29 S$=0-047 C8=2=N.613 LYRB=6.86 B/T=2-81 FULLY LOADED 
15.5 0.047 e7.0 27.0 0.613 0.0 12.0 
0,0 
0.90280 0.095 0.019] 0.069 
2 MARINER TIME (SEC) FOQCE (KIPS) 
AL ANK 
BLANK 


4/7 END 


FIGURE 26a - 


SAMPLE PROGRAM - DATA DECK 
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THIS PROGRAM CALCULATES THE SLAM FORCE DISTRIBUTION vS TIME 
ON THE LOWER 1/18 OF THE SRIP FROM STATION 8 FORWARD. 


MARTNER © = w= = we © © @ @ ow © ow oo @ AS BUILT 

WALF BEAM 38.0aCAa DFSIGN DRAFT s 29¢75000 

, ee = 5PR.AD LENSTH LOAD WATER LINE s 537-13 
DISTANCE BETWEFN STATIONS = Pane 


(ALL DISTANCES IN FEET) 


SSSTSH SS eSSE SSH HCAS ESET SCE SS ESE eEeeEaneOeaseneee 


» es 
a AVERAGE HULL SMAFE FACTOR » 9011909 s 
e s 


SOSH SK SCH HHCSHSH Ke TEKST SHES SE SEE EHH EE HORE EEE 


FIGURE 26b - SAMPLE PROGRAM_- OUTPUT 
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